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Abstract 

A detailed study of pseudorapidity densities and multiplicity distributions of primary charged parti¬ 
cles produced in proton-proton collisions, at ^fs = 0.9, 2.36, 2.76, 7 and 8 TeV, in the pseudorapidity 
range I;;! < 2, was carried out using the ALICE detector. Measurements were obtained for three event 
classes; inelastic, non-single diffractive and events with at least one charged particle in the pseudo¬ 
rapidity interval It;! < 1. The use of an improved track-counting algorithm combined with ALICE’S 
measurements of diffractive processes allows a higher precision compared to our previous publica¬ 
tions. A KNO scaling study was performed in the pseudorapidity intervals I 77 I < 0.5, 1.0 and 1.5. The 
data are compared to other experimental results and to models as implemented in Monte Carlo event 
generators PHOJET and recent tunes of PYTH1A6, PYTHIAS and EPOS. 
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1 Introduction 

The multiplicity of emitted charged particles is one of the most basic characteristics of high-energy 
hadron collisions and has been the subject of longstanding experimental and theoretical studies, which 
have shaped the understanding of the strong interaction. Following on from earlier ALICE studies of 
global properties of proton-proton (pp) collisions l[T]® , this publication presents a comprehensive set of 
measurements of the pseudorapidity density (dAch/d?/) of primarjQ charged particles and of their multi¬ 
plicity distributions over the energy range covered by the LHC, from 0.9 to 8 TeV. The pseudorapidity 
density of primary charged particles was studied over the pseudorapidity range \r]\ < 2 , and their multipli¬ 
city distributions in three intervals: \r]\ < 0.5, 1.0 and 1.5. Results are given for three conventional event 
classes: (a) inelastic (INEL) events, (b) non-single diffractive (NSD) events and (c) events with at least 
one charged particle in I 77 I < 1 (INEL>0). 

At LHC energies, particle production is still dominated by soft processes but receives significant contri¬ 
butions from hard scattering, thus multiplicity and other global event properties measurements allow to 
explore both components. As these properties are used as input in Glauber inspired models ll9UT2]l. such 
studies are also contributing to abetter modelling of Pb-Pb collisions. Already at 8 TeV, high multiplicity 
proton-proton collisions provide energy densities comparable, for instance, to energy densities in Au-Au 
central collisions at RHIC, allowing a comparison of nuclear matter properties in strongly interacting 
systems with similar energy densities but with volumes orders of magnitude smaller. 

It is worth noting that, already - 2.36 TeV, hadron collision models tuned to pre-LHC data failed 
to reproduce basic characteristics of proton-proton collisions at the LHC, such as pseudorapidity den¬ 
sity of charged particles, multiplicity distributions, particle composition, strangeness content, transverse 
momentum distributions and sphericity (see for instance [2], Q, ||11 and 113). Therefore, a more pre¬ 
cise measurement of charged-particle multiplicity distributions and a study of their energy dependence 
contribute to a better understanding of particle production mechanisms and serve to improve models. In 
turn, a better simulation of collision properties improves the determination of the detector response and 
background estimates of underlying event properties relevant to the study of high-px phenomena. 

In the Regge theory 1141 . one of the most successful models for describing soft hadronic interactions, 
the asymptotic behaviour of cross-sections for elastic scattering and multiple production of hadrons is 
determined by the properties of the Pomeron, the t-channel right-most pole, in the elastic scattering 
amplitude. In QCD, the Pomeron, which has vacuum quantum numbers, is usually related to gluonic 
exchanges in the t-channel. The experimentally observed increase of the total cross-section with increas¬ 
ing collision energy made it necessary to consider a Pomeron as a Regge trajectory with t = 0 intercept: 
a/)(0) = 1 -I- A > 1 1141 . The energy dependence of the particle (pseudo-)rapidity density provides in¬ 
formation about the Pomeron trajectory intercept parameter, A. If interactions between Pomerons are 
neglected, the inclusive particle production cross-section, crinci., is determined only by the contribution 
of the single (cut-)Pomeron exchange diagram. In this approximation, dcrind./dy (~ dcrinci./dp) at mid¬ 
rapidity is proportional to lITSll . Thus, the energy dependence of the inclusive cross-section gives 
more reliable information about the value of A than the energy dependence of the total interaction cross- 
section, for which contributions from multi-Pomeron exchanges strongly modify the energy dependence 
of the single Pomeron exchange diagram. In the same approximation, the energy dependence of the 
particle (pseudo-)rapidity density in the central rapidity region is given by dV/dy oc s^j crim., where cint. 
is the interaction cross-section (see for instance llT6l - IT9]l l. Up to LHC energies, crim. is well represented 
by a power law of s. However, for reasons of unitarity EOll . it is expected that this power law should be 
broken at sufficiently high energy, although well above LHC energies. Therefore, the energy dependence 
of the particle (pseudo-) rapidity density in the central region at LHC, dV/dy w dNjAr], should follow 
the same power law trend. In this publication, this relationship is explored further for three event classes 

^Primary particles are defined as prompt particles produced in the collision including all decay products, except those from 
weak decays of light flavour hadrons and muons. 
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and using 5 ALICE data points. 

It was more than 40 years ago that A. M. Polyakov ll^ and then Z. Koba, H.B. Nielsen and P. Olesen Il22l 
proposed that the probability distribution of producing n particles in a collision, P{n), when expressed as 
a function of the average multiplicity, (n), should reach an asymptotic shape at sufficiently high energy 



( 1 ) 


where 'P is a function supposed to describe the energy-invariant shape of the multiplicity distribution. 


Such scaling behaviour is a property of particle multiplicity distributions known today as Koba-Nielsen- 
Olesen (KNO) scaling. 

One well identified mechanism for KNO scaling violafion is fhe increasing probability of multi-parton 
scattering with increasing ^/s. Moreover, since the topologies and multiplicities of diffractive and non- 
diffractive (ND) events are different, their KNO behavior may be different. Even if KNO scaling were 
to be valid for each, it might not be valid for their sum. Nevertheless, KNO scaling is expected to be 
violated for both diffractive and non-diffractive processes l24l at sufficiently high collision energies 
and the LHC provides the best opportunity to study the extent of these scaling violations. 

Indeed, deviation from KNO scaling was already observed long ago at ISR energies (proton-proton 
collisions atV^ from 30.4 to 62.2 GeV), in the full phase space, for inelastic events ESi . On the other 
hand, for NSD collisions, scaling was still found to be present E5l . suggesting that diffractive processes 
might also play a role in KNO scaling violations. In e'''e“ collisions, at from 5 to 34 GeV, KNO scaling 
was found to hold within ±20% E^ . In proton-antiproton collisions at the CERN collider = 200, 
546 and 900 GeV), KNO scaling was found to be violated for NSD collisions in full phase space ETll . 
ESll . E9ll . Nevertheless, for NSD collisions, in limited central pseudorapidity intervals, KNO scaling 
was still found to hold up to 900 GeV, and at = 546 GeV, KNO scaling was found to hold in the 
pseudorapidity interval \r]\ < 3.5 BUIETI . In NSD proton-proton collisions at the LHC, at- 2.36 and 
7 TeV and in I 77 I < 0.5, ALICE f2l and CMS E2]| observed no significant deviation from KNO scaling. 

This publication presents a study of KNO scaling, at from 0.9 to 8 TeV, in three pseudorapidity inter¬ 
vals (|? 7 | < 0.5, 1.0 and 1.5) and for a higher multiplicity reach compared to previous ALICE publications, 
quantified wifh KNO variables (momenfs) E2l as well as wifh fhe paramefers of Negative Binomial Dis- 
fribufions (NBD) used fo fif measured mulfiplicify disfribufions. 

Wifh respecf fo previous ALICE publicafions, fhe analysis reported here makes use of improved fracking 
and frack-counfing algorifhms; beffer knowledge and improved simulafion of diffracfion processes; an 
expanded pseudorapidify range for dVch/d/y sfudies and beffer sfafisfical precision af - 0.9 and 7 TeV, 
extending by a factor of 2 the previously published multiplicity distribution reach. Results at = 2.76 
and 8 TeV are presented for the first time in this publication. 

Previous measurements of both dVch/d?/ and multiplicity distributions from CMS E3ll34l and UA5 ETl 
allow a direct comparison to our data. Others by ATLAS 051 and LHCb E6]l use different definitions (77 
and pt ranges) making direct comparison impossible. 

This publication is organized as follows: Section describes the ALICE sub-detectors relevant to this 
study; Section provides the details of the experimental conditions and of the collection of data; Sec¬ 
tion 1^ explains the event selection; Section describes the track selection criteria and the three track 
counting algorithms; Sections and report the analyses for the measurement of the pseudorapidity 
density and of multiplicity distributions, respectively; Section [^discusses systematic uncertainties; Sec¬ 
tion [^presents the multiplicity measurements, NBD fits of the multiplicity distributions, KNO scaling 
and ^-moment studies. Linally, in Sectionthe results are summarized and conclusions are given. 
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2 ALICE subdetectors 

The ALICE detector is fully described in |[37ll . Only the main properties of subdetectors used in this 
analysis are summarized here. Charged-particle tracking and momentum measurement are based on 
data recorded with the Inner Tracking System (ITS) combined with the Time Projection Chamber (TPC) 
(381, all located in the central barrel of the ALICE detector and operated inside a large solenoid magnet 
providing a uniform 0.5 T magnetic field parallel to the beam line. 

The VO detector (Wl consists of two scintillator hodoscopes, each one placed at either side of the inter¬ 
action region, at z = 3.3 m (VOA) and at z = -0.9 m (VOC) (z is the coordinate along the beam line, with 
its origin at the centre of the ALICE barrel detectors), covering the pseudorapidity ranges 2.8 < 77 < 5.1 
and -3.7 < 77 < -1.7, respectively. The time resolution of each hodoscope is better than 0.5 ns. 

The ITS is composed of high resolution silicon tracking detectors, arranged in six cylindrical layers at 
radial distances to the beam line from 3.9 to 43 cm. Three different technologies are employed. Eor 
the two innermost layers, silicon pixels (SPD (401 ) are used, covering pseudorapidity ranges I 77 I < 2 and 
I 77 I < 1.4, respectively. The SPD is followed by two Silicon Drift Detector layers (SDD, iHdl i. The Silicon 
Strip Detector (SSD, (421 ) constitutes the two outmost layers consisting of double-sided silicon micro¬ 
strip sensors. The intrinsic spatial resolution (cr^^ x cr^) of the ITS subdetectors is: 12x100 /im^ for SPD, 
35x25 for SDD, and 20x830 for SSD, where tp is the azimuthal angle and r the distance to 
the beam line. The ITS sensors were aligned using survey measurements, cosmic muons and collision 
data (431 . The estimated alignment accuracy is 8 77 m for SPD and 15 /rm for SSD in the most precise 
coordinate {np). Eor the SDD, the intrinsic space point resolution is cr^ = 30 yum in the z direction and 
o'rtf = 40 to 60 77 m, depending on the sensor, along rp (drift). Because of some anomalous drift field 
disfribufions, in fhe reconsfrucfion, a sysfemafic uncerfainfy up fo 50 77 m in z and 500 77 m in rp was 
added fo accounf for differences befween dafa and simulafion. The ITS resolufion in fhe deferminafion 
of fhe fransverse impacf paramefer measured wifh respecf fo fhe primary verfex is fypically 70 77 m for 
Tracks wifh px = 1 GeV/c, including fhe confribufion from fhe primary verfex posifion resolufion. 

The SPD and fhe VO scinfillafor hodoscopes provided Triggers for collecting dafa. 

The TPC (38l is a large cylindrical driff defecfor wifh a cenfral high volfage membrane af z = 0, main- 
fained af -i-lOO kV and Two readouf planes af fhe end-caps. The maferial budgef befween fhe inferacfion 
poinf and fhe active volume of fhe TPC corresponds fo 11% of a radiafion lengfh, when averaged over 
|77|<0.8. 

The TPC and fhe ITS were aligned relafive fo each ofher wifhin a few hundred micromefers using cosmic- 
ray and proton collision dafa (431 . 

The momenfum measuremenf is nof explicifly used in This sfudy, however, fhe simulafion of fhe defecfor 
response is sensifive fo fhe particle momenfum specfrum. Since evenf generafors used in Monfe Carlo 
simulations do nof reproduce fhe observed momenfum disfribufions, fhe difference befween dafa and 
Monfe Carlo simulafion is Taken info accounf when evaluating sysfemafic errors. Eor momenfa lower Than 
2 GeV/c, represenfing fhe bulk of fhe dafa, fhe 777 resolufion for Tracks measured in fhe TPC and in fhe 
ITS, is abouf 0.80% af 777 = 1 GeV/c, if increases fo 0.85% af 777 = 2 GeV/c and fo 3% af 777 = 0.1 GeV/c. 

Charged-parficle mulfiplicifies were measured using informafion from fhe TPC in I 77 I < 0.9 and from fhe 
ITS in I 77 I < 1.3. Af larger pseudorapidifies, fhe SPD alone was used to expand fhe range of dVch/d 77 
measuremenfs fo I 77 I < 2 . 0 . 
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3 Experimental conditions and data collection 

3.1 Proton beam characteristics 

Data were selected during LHC collision periods at a luminosity low enough to allow the minimum bias 
trigger rate not to exceed 1 kHz. At = 0.9 TeV, the number of protons per colliding bunch varied 
from 9x10^ to 3.4xl0'\ while the number of colliding bunches was either 1 or 8. At = 2.76 TeV, 
the number of protons per colliding bunch varied from 5x10'^ to 7x10^^, while the number of colliding 
bunches was either 48 or 64. At = 7 TeV, the number of protons per colliding bunch varied from 
8.6x10^ to 1.4x10^^, resulting in a luminosity between 10^^ and 10^° cm“^s“'. There were up to 36 
bunches per beam colliding at the ALICE interaction point. When needed, the luminosity was kept below 
lO^** cm“^s“^ by a transverse displacement of the beams with respect to one another. AtV^ - 8 TeV, there 
were 3 proton bunches colliding at the ALICE interaction point each containing about 1.6x10^^ protons. 

Data used for this study were collected at low beam currents, so that beam-induced backgrounds (beam- 
gas or beam-halo events) were low and could be removed offline using VO and SPD detector information, 
as discussed in Section l4~n 

3.2 Triggers 

The ALICE trigger system is described in 14411 . Data were collected with a minimum bias trigger, MBqr, 
requiring a hit in the SPD or in either one of the VO hodoscopes; i.e. essentially at least one charged 
particle anywhere in the 8 units of pseudorapidity covered by these detectors. Triggers were required 
to be in time coincidence with a bunch crossing the ALICE interaction point. Control triggers, taken 
for various combinations of beam and empty-beam buckets, were used to measure beam-induced and 
accidental backgrounds. 

3.3 Characteristics of data samples used in this study 

General characteristics of the data samples used are given in Table [T] 

The data at = 0.9 TeV were collected in May 2010, with one polarity of the ALICE solenoid magnet 
(solenoid magnet field pointing in the positive z direction). 

The first LHC data above Tevatron energy were collected in 2009, at = 2.36 TeV, in a run with unstable 
LHC beams, during which only the SPD was turned on. Therefore, in this case, the charged-particle 
multiplicity was measured using exclusively the SPD information. In this publication, the previously 
published results at = 2.36 TeV 0 are used for comparison. 

Proton-proton data were collected at = 2.76 TeV, an energy that matches the nucleon-nucleon centre- 
of-mass energy in the first Pb-Pb collisions provided by the LHC, in 2011. 

Data at = 7 TeV were collected in 2010. About 20% of the data were taken with a magnet polarity 
opposite (solenoid field poinfing in fhe negative z direction) fo fhaf of^/s = 0.9 TeV dafa. A sample of 
12.3x10^ evenfs, collecfed wifhouf magnetic field, was used fo check some of fhe sysfemafic biases in 
frack reconsfrucfion. 

Af = 8 TeV only a subsef of runs was collecfed wifh fhe MBqr as a minimum bias frigger in 2012, 10 
were selecfed for fhis analysis. 

Af 0.9 and 7 TeV, dafa samples are subsfanlially larger fhan fhose available in previous ALICE publi- 
cafions on charged-parficle mulfiplicifies UHll. Lor fhe charged-parficle mulfiplicify analysis, fhe evenf 
sample af = 0.9 and 7 TeV increased by a factor of 50 and 2000, respectively, giving significanf 
exfension of fhe mulfiplicify reach and beffer sfafisfical precision. The precision of dVch/d ?7 is nof sub- 
sfanfially limifed by evenf sample size. However, fhe large number of runs available made if possible fo 
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MB events (xlO^) 


Vs (TeV) 

Triggered 

Reconstructed 

Selected 

(B) 

Luminosity (nb V 

0.9 

7.4 

6.3 

5.6 

0.04 ±0.01 

0.128 ±0.006 

2.36 

0.097 

0.097 

0.04 

< 0.001 

- 

2.76 

33.9 

32.6 

28.3 

0.025 ±0.01 

0.583 ±0.013 

7 

404.4 

384.2 

343.7 

0.04 ±0.01 

6.05 ±0.25 

8 

31.5 

26.6 

24.1 

0.02 ± 0.01 

0.41 ±0.02 


Table 1: For each centre-of-mass energy: total number of minimum bias (MB) events collected; number of those 
events that were reconstructed; number of reconstructed events passing the selection described in the text, except 
for z vertex quality and position; average number of interactions per bunch crossing, (//); integrated luminosity 
corresponding to the number of events reconstructed. 

study run-to-run fluctuations of the dAch/d?? measurements over long periods of time, thus providing a 
monitoring of the uniformity of the data quality. 

4 Event selection 

4.1 Background rejection 

4.1.1 Beam background 

The main sources of event background are beam gas and beam halo collisions. Such events were removed 
by requiring that the timing signals from the VO hodoscopes, if present, be compatible with the arrival 
time of particles produced in collision events. In addition, because of the different topology of beam 
background events, the ratio between the number of SPD clusters and the number of SPD tracklet^is 
much higher in beam background events, therefore a cut on this ratio was applied. The remaining fraction 
of beam background events in the data, estimated by analysing special triggers taken with non-colliding 
bunches or empty beam buckets, does not exceed 10“'^ for all centre-of-mass energies. The track beam 
background is mostly significant in the last 77 bins (I 77 I w 2 ) where it reaches 4x10“^ in the worst case. 

4.1.2 Event pileup 

The other type of potential event background comes from multiple collision overlap. For the data 
used in this publication, the proton bunch spacing was 50 ns or longer, the luminosity did not exceed 
10^** cm“^s“', and the probability to have collisions from different bunch crossings in the 300 ns inte¬ 
gration time of the SPD was negligible. However, multiple collisions in the same bunch crossing, also 
referred to as event pileup or overlap, have to be considered in case their vertices are not distinguish¬ 
able. In order to avoid or minimize corrections for event pileup, runs with a low number of interactions 
per bunch crossing, p < 0.061, were selected resulting in an average p, (p) < 0.04, for all data samples 
(Table [T]). This corresponds to at most 2% probability of more than one interaction per event. 

The identification of pileup events relies on multiple vertex reconstruction in the SPD, with algorithms 
using three basic parameters: (a) The distance of closest approach (DCA) to the main vertex for a SPD 
tracklet to be included in the search for an additional interaction: DCA > 1 mm; (b) The distance between 
an additional vertex and the main vertex, Az > 8 mm; (c) The number of SPD tracklets (Atrk) used to 
determine an additional vertex (number of contributors to the vertex): Atrk > 3. 

With this choice of parameters, and with the relatively broad z vertex distribution at the LHC (FWHM > 
12 cm), typically only 10% to 15% of multiple collisions are missed, and the fraction of fake multiple 
collisions due to SPD vertex splitting from a single interaction is low (typically a few times 10~^). 

tracklet is a short track segments in the SPD, compatible with the event vertex. 
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The pileup detection efficiency was studied both by overlapping two Monte Carlo proton-proton colli¬ 
sions and by measuring pileup in the data. The pileup fraction, estimated from identified pileup events 
in the data, is found to be consistent with what is expected from the values derived from trigger infor¬ 
mation (Table [T]l. 

In multiplicity measurements, pileup affects the data mainly when two vertices are not distinguishable. 
When they are distinguishable, the multiplicity is taken from the vertex with the highest number of tracks. 
The small bias induced by choosing systematically the highest multiplicity vertex is negligible in our low 
pileup data samples. 

Comparing dAch/d?/ measurements, for different runs, no correlation is found between dAch/d?/ values at 
77 = 0 and /I values. Comparing data with and without identified pileup rejection, the change in dAch/d?/ 
values is smaller than 0.5%, which is smaller than systematic uncertainties. Note that the requirements 
for track association to the main vertex reject a further fraction of the tracks coming from the 10% to 
15% of unidentified pileup collisions. The conclusion is that event pileup corrections to dAch/d?/ are 
negligible in these low pileup data samples. 

For multiplicity distributions, even though data were selected with a low pileup probability, it is impor¬ 
tant to verify that the pileup does not distort the distributions, as the relative pileup fraction increases 
with multiplicity. The fraction of pileup events, which the ALICE pileup detection algorithm identifies 
after the event selection, is about 10“^, with no significant differences between the four centre-of-mass 
energies. Moreover, tight DCA cuts allow tracks originating from the main vertex to be distinguished 
from those coming from a pileup vertex even when the vertices are closer than 0.8 cm in z. This was 
confirmed by simulating events, where two Monte Carlo pp collisions were superimposed, demonstrat¬ 
ing that only 5% of the events passing the selection had extra tracks from the secondary vertex. In 90% 
of such cases, the distance along the beam line between the two vertices was Az < 0.5 cm. In the data 
samples with a pileup fraction of order yu/2 < 0.02, the residual average fractions of events with pileup 
is at most 0.4%. Furthermore, the simulation shows that the pileup that does affect the multiplicity of an 
event is rather broadly distributed across events with different multiplicity, but becomes significant only 
outside the multiplicity range studied here. The multiplicity at which the pileup contribution reaches 
10% of the measured multiplicity at = 7 TeV is Ach = 105, 170 and 310, for \r]\ < 0.5, 1.0, and 1.5, 
respectively, which is beyond multiplicity ranges covered in this publication. 

Therefore, no pileup corrections were applied. Other background contributions from cosmic muons or 
electronics noise are also negligible. 

4.2 Offline trigger requirement 

Both for the INEL and INEL>0 normalizations, the online MBqr trigger was used. However, for the 
NSD analysis, a subset of the total sample was selected offline by requiring a coincidence (MB and) 
between the two VO hodoscope arrays. This corresponds to the detection of at least one charged particle 
in both hemispheres, in the VO hodoscope arrays separated by 4.5 units of pseudorapidity, a topology that 
tends to suppress single-diffraction (SD) events; therefore, model dependent corrections and associated 
systematic errors are minimized. 


4.3 Vertex requirement 


The position of the interaction vertex is obtained either by correlating hits in the two silicon-pixel layers 
(SPD vertex), or from the distribution of the impact parameters of reconstructed global track^ (global 
track vertex) llT7l[^l4^l45ll4^ . The next step in the event selection consists of requiring the existence 
of a reconstructed vertex. 


^Tracks reconstructed in the TPC and matched to ITS clusters (see Section 5.1 and references therein) 
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Fig. 1 : For a data sample of events fulfilling the MBqr trigger selection, at Vs = 7 TeV, the distribution of the 
quantity d^Nchj dz/dz is plotted for tracklets, in the plane pseudorapidity ( 77 ) vs. z position of the SPD vertex (Zvtx). 
showing the dependence of the 77 acceptance on Zvtx- 


Two SPD vertex algorithms were used: a three-dimensional vertexer (3D-vertexer) that reconstructs the 
X, y and z positions of the vertex, or a one-dimensional vertexer (ID-vertexer) that reconstructs the z 
position of the vertex. The vertex position resolution achieved depends on the track multiplicity. For the 
3D-vertexer it is typically 0.3 mm both in the longitudinal (z) direction and in the plane perpendicular to 
the beam direction. The ID-vertexer resolution in the z direction is on average 30 /urn. If the 3D-vertexer 
algorithm does not find a vertex (typically 47% of the cases at V^ = 7 TeV), then the simpler ID-vertexer 
is used to determine the z position of the vertex, and the x and y coordinates are taken from the average 
X and y vertex positions of the run. The 3D-vertexer efficiency is strongly multiplicity dependent. As 
the bulk of the events have a low multiplicity, this explains the relatively low average vertex finding 
efficiency. For fhe z coordinafe, if no reliable verfex is found (fypically 14% of fhe cases), eifher because 
the ID-vertexer did not find a verfex or fhe ID-verfex qualify was nof sufficienf (fhe dispersion of fhe 
difference of azimufhal angles befween fhe fwo hifs, one in each SPD layer, of fracklefs confribufing fo 
fhe verfex is required fo be smaller fhan 0.02 rad), fhe evenf is rejecfed. For fhe global frack verfex, fhe 
resolution is fypically 0.1 mm in fhe longifudinal (z) direction and 0.05 mm in fhe direcfion fransverse fo 
fhe beam line. 

Bofh SPD and global frack vertices have fo be presenf and consisfenf by requiring fhaf fhe difference 
befween fhe fwo z positions be smaller fhan 0.5 cm. If nof, in 3 fo 4% of fhe cases, fhe evenf is rejecfed. 
The cuf was chosen fo be compafible wifh DCAj, cuf applied fo fracks fo ensure fhaf we combine fracklefs 
and fracks from fhe same collision (see Secfio n|5.2| ). This condifion removes mainly non-Gaussian fails 
in fhe columns of fhe defecfor response mafri)(|^ af low mulfiplicify, coming from fhe facf fhaf SPD and 
frack vertices, when separafed, fend fo have differenf mulfiplicifies associafed fo fhem. In fhe dafa. Ibis 
requiremenf also removes 80% of pileup evenfs wifh well-separafed vertices. 

The T] accepfance is correlafed wifh fhe verfex z posifion (Zvtx) (Fig. [^. For mulfiplicify disfribufion 
measuremenfs, in order for fracks fo remain wifhin fhe accepfance of fhe SPD in fhe 77 versus Zvtx plane, 
fhe following requiremenfs were imposed on fhe vertex posifion along fhe z axis: Izvtxl < 10, 5.5 and 1.5 
cm for |? 7 | < 0.5, 1 and 1.5, respecfively. In fhe measuremenf of dVch/d??, fhe requiremenf on fhe verfex 
was relaxed fo |zvtxl < 30 cm, in order fo allow exfending fhe 77 range fo I 77 I < 2. 

^The response matrix is a 2-dimensional matrix obtained from the simulation giving the correspondence between generated 
and observed multiplicities. A response matrix column consists of the digitized distribution of the probability to measure a 
given multiplicity for a given generated multiplicity. 
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4.4 Event selection efficiency 

As described in BTI . PYTHIA6 HStiSOl and PHOJET 115111521 event generators used by AEICE were 
adjusted to reproduce the measured diffraction cross-sections and the shapes of the diffracted mass {Mx) 
distributions extracted from the Kaidalov-Poghosyan model |[5Jl . These modified versions of evenf gen- 
erafors are referred fo as “funed for diffracfion”. Typically, ctsd/cinel ~ 0.20, where ctinel is the inelas¬ 
tic cross-section, ctsd is the SD cross-section for Mx < 200 GeV/c^, and ctdd/ctinel » 0.11, where ctdd 
is the double diffraction cross-section for A 77 > 3 (A?/ is the size of the particle gap in the pseudorapidity 
distribution). These fractions have insignificant energy dependence between 0.9 and 7 TeV 1471 . and the 
values at 7 TeV were used for 8 TeV data. 


Table [T] shows the number of events selected at each centre-of-mass energy prior to the Zytx requirement. 
Selection efficiencies using criferia defined above in fhis section, were esfimafed for INEE, NSD and 
SD evenfs (classified af generafor level by evenf generafor flags) as a funcfion of fhe number of gen- 
erafed charged parficles (shown on Eig. for the case \r]\ < 1 and the various centre-of-mass energies 
considered). The particular selection is designated by the offline trigger used to construct it, MBqr or 
MB and- Note that for dVch/d ?7 measurement selection efficiencies are defined in a separafe way (see 
Secfion 6.11. TeV fhe INEE even! selecfion efficiency based on fhe MBqr trigger reaches 

100 % for a charged-particle multiplicity above 8. 


Eor SD events, the efficiency of the MB and selection reduces significantly when going to higher energies 
(Fig.ig, because the Eorentz boost of the diffracted system increases with increasing centre-of-mass 
energies. This implies that in the normalization to the NSD event class, corrections for the remaining SD 
contribution become smaller when going to higher energies. The MB and trigger selects 84%, 86%, 87% 
and 87% of the MBqr triggers, and 13%, 4%, 1% and 1 % of the SD events satisfy the MB and selection, 
atV^ = 0.9, 2.76, 7 and 8 TeV, respectively. 


5 Track selection and multiplicity algorithms 
5.1 Track quality requirements 

The following criteria were used to select reconstructed tracks associated to the main event vertex: 


- for tracks reconstructed from both ITS and TPC information (global tracks), the selection requires at 
least 70 pad hit clusters in the TPC, a good track quality dof < 4), a distance of closest approach 
(DCA) along the z direction (DCA^) < 0.5 cm, and a pT-dependent transverse DCA (DCAt) require¬ 
ment, which corresponds to a 7 sigma selection. DCAt conditions are relaxed by a factor 1.5 for tracks 
lacking SPD hits. 


- for tracks reconstructed with ITS information only (ITS-only tracks) the number of ITS hit clusters 
associated to the track must be larger than 3, among the 6 layers of the ITS, and x^j dof < 2.5. The 
DCA^ and DCAt requirements are the same as for global tracks. 

- for SPD tracklets, the association to the vertex is ensured through &x^ requirement. Using the SPD 
vertex as the origin, differences in azimuthal (A(^ = (p 2 -<P\, bending plane) and polar (A0 -92-0\, non¬ 
bending direction) angles are calculated between hits in the inner (layer 1) and in the outer (layer 2) 
SPD layers. Hit combinations, called tracklets, are selected with the following condition 






CTTr 


2 f Si +02 \ 


sm 


{Aey 


err, 


< 1.6 


( 2 ) 


where = 0.08 rad, ere - 0.025 rad and the sin^ factor takes into account the 6 dependence of A9. 
The x^ value 1.6 was chosen to lie well within the part of the x^ distribution of the data correctly 
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Fig. 2: Charged-particle multiplicity (Nch - number of primary charged particles generated in I 77 I < 1) dependence 
of the efficiency of the event selection described in Section]^ obtained as the average between PYTHIA6 PerugiaO 
and PHOJET, both tuned for single diffraction defined for Mx < 200 GeV/c^ (see Il47l i. Efficiencies are given for 
INEL events with MBqr trigger (open circles), NSD events with MB and trigger (open squares), and SD events 
with MB and trigger (open diamonds), — 0.9 TeV (top left), 2.76 TeV (top right), 7 TeV (bottom left) and 
8 TeV (bottom right). Error bars correspond to the difference between the two event generators and statistical 
uncertainty added in quadrature (non-negligible only for the SD events selection efficiency). 

reproduced by the simulation. The cut imposed on the difference in azimuthal angles rejects charged 
particles with a transverse momentum below 30 MeV/c; however, the effective transverse-momentum 
cut-off is defermined mosfly by parficle absorpfion in fhe maferial and is approximafely 50 MeV/c, in 
|? 7 l < 1. If more fhan one hif in an SPD layer mafches a hif in fhe ofher layer, only fhe hif combination 
wifh fhe smallest value is used. 

Some of the SPD elements had to be turned off, resulfing in lower efficiency in some regions of fhe 
77 versus azimufhal angle plane. In order fo reach fhe besf possible precision in fhe measuremenf of 
dAch/d/y, fiducial cufs were applied fo bofh fracks and fracklets, excluding azimufhal regions where fhe 
fracking efficiency correcfions are relatively large. These fiducial cufs vary wifh dafa faking periods, 
following fhe evolution of fhe SPD accepfance. Af = 0.9, 2.76, 7, and 8 TeV, fhe fractions of fhe 
accepfance removed were 64%, 68 %, 65%, and 35%, respectively. Some of the SPD elements could be 
recovered before collecting 8 TeV data, explaining the improvement. 

For multiplicity distribution studies, fiducial cuts were not applied because they increase statistical un¬ 
certainty, hence limiting the high multiplicity reach. 

5.2 Track counting algorithms 

In previous ALICE publications |[T]j3l> the charged-particle multiplicity was measured in |? 7 | <1.3 using 
only SPD tracklets built from SPD pixel hits. In order to extend the pseudorapidity range to I 77 I < 2, an 
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improved tracklet algorithm, initially used in (541, was introduced to take into account the 9 dependence 
of the uncertainty in thex^ (Eq. ([^). With this improvement, the efficiency for detecting SPD tracklets 
became uniform as a function of pseudorapidity and z position of the vertex, which allowed vertices 
further away from the nominal interaction point along the beam direction to be used, thereby extending 
significantly the pseudorapidity range. 

To be less sensitive to the SPD acceptance, track counting algorithms were developed, that make use of 
tracking information from other AEICE detectors, the SDD, the SSD and the TPC. Each track is counted 
as primary if it fulfills fhe transverse DC A requirements listed in Section [5A] and it is not associated to a 
secondary vertex identified by a dedicafed algorifhm (451 funed fo fag y-conversions, and A decays. 

Three mulfiplicify esfimafors were developed by AEICE using fhree differenf samples of fracks: 


- SPD fracklefs, wifh \'q\ < 2 (referred fo as Tracklef algorifhm^ The Tracklef algorifhm sfores, for each 
fracklef, references fo ITS or global frack candidafes using af leasf one of ifs pixel clusfers. 

- ITS-only fracks, wifh \t]\ < 1.3, obfained using all hil clusfers in fhis defecfor, plus fracklefs (|? 7 | < 2) 
built out of SPD pixel clusters not matched to any ITS track (referred to as ITS+ algorithm). 

- TPC tracks, with \r]\ < 0.9, matched to hits in the ITS, plus ITS-only tracks (up to \ri\ < 1.3) built out of 
silicon hit clusters not matched to any TPC track, plus tracklets (It/I < 2) built out of SPD pixel clusters 
not matched to any ITS or TPC track (referred to as ITSTPC-t algorithm). 
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Fig. 3: Graphical representation of the detector response matrices obtained with PYTHIA6 CSC (481 combined 
with a simulation of the ALICE detector, at = 7 TeV, for three pseudorapidity intervals (It/I < 0.5, 1.0, and 1.5 
from left to right, respectively), and for the three track counting algorithms, Tracklet, ITS-h and ITSTPC-r, from 
top to bottom, respectively. Horizontal axes show generated primary charged-particle multiplicities and vertical 
axes measured multiplicities. 


In order to keep away from the edges of the detectors, where the acceptance is less precisely known, ITS 
and TPC tracks used in this study are limited to |t/| <1.3 and |t/| < 0.9, respectively. 

^Potentially |;;| < 3 can be reached using event vertices displaced from the detector center at distances Izvtxl ^ 30 cm (see 
Fig. [Tf, however the sample size of such events is too small. 
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Fig. 4: The three multiplicity algorithms are compared, after full correction, (left) without and (right) with fiducial 
cuts in azimuthal angles. Ratios of dNch/di] measurements with different algorithms are shown: ITSTPC+ over 
Tracklet (black circles) and ITS+ over Tracklet (red squares). 

Properties of the three track counting algorithms are compared in Fig. showing that, going from 
Tracklet to ITS+ and to ITSTPC+ algorithms, the detector response matrix becomes narrower and has 
a topology closer to that of a diagonal matrix. When going from \r]\ < 0.5 to \t]\ < 1.5, the response 
matrix becomes broader and has a less diagonal topology, as geometrical acceptance effects become 
more important, and dominated by the SPD with significant inefficiency due fo some missing modules. 
Nofe fhaf by resfricfing fhe azimufhal angle fo good regions of fhe SPD, fhe difference befween algorifhms 
in dAch/d?/ measuremenfs is of order ±1% in fhe cenfral region (Fig. |^. However, fhe resulf wifh fhe 
Tracklef algorifhm is nof sensifive fo fhis cuf, and, as if is needed fo measure mulfiplicifies beyond 
\t]\ = 1.3, if is used alone for dAch/d?/ measuremenf. For mulfiplicify disfribufion measuremenfs all 
fhree algorifhms are used wifhouf fhe ip region resfricfions wifh a corresponding sysfemafic uncerfainfy 
confribufion. 

In fhe pseudorapidify region I 77 I < 0.9, fhe TPC accounfs for 90% of fhe fracks, fhe ITS complemenf 9% 
and fhe SPD complemenf 1%. These fractions vary wifh fhe rj range. Oufside I 77 I < 1.3, SPD fracklefs are 
fhe only confribufion. The small flucfuafions befween poinfs in \r]\ >1.3 come from fhe slighfly differenf 
number of evenfs used for averaging befween algorifhms, affer efficiency correcfions in each 77 bin. 

6 Pseudorapidity density of primary charged particles: analysis 

Raw dAch/d?/ disfribufions have fo be correcfed for defecfor and trigger acceptance and efficiency, and 
for confaminafion from daughfers of sfrange particles. Nofe fhaf fhis secfion only describes fhe parficu- 
larifies of dAch/d?/ measuremenf, unless specifically sfafed ofherwise. For charged particle mulfiplicify 
disfribufion measuremenf see Secfion |7] 

6.1 Acceptance and efficiency corrections 

Three fypes of correcfions have fo be applied fo fhe raw dafa: (a) a frack-fo-parficle correcfion fo fake info 
accounf fhe difference befween measured fracks and “frue” charged primary particles. This correcfion 
mainly depends on accepfance effecfs and on defecfor and reconsfrucfion efficiencies; (b) correcfions 
for fhe bias coming from fhe vertex reconsfrucfion requiremenf, af bofh frack and evenf levels (verfex 
reconsfrucfion correcfion). This bias exisfs on bofh fhe number of fracks and fhe evenfs used, since 
evenfs wifhouf a reconsfrucfed verfex are nof selecfed, and fracks from fhose evenfs fherefore do nof 
confribufe; (c) correcfions af bofh frack and evenf levels, fo fake info accounf fhe bias due fo fhe MBqr 
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trigger required for INEL and INEE>0 event classes or the MB and offline selection for the NSD event 
class. 

In practice, the number of tracks is corrected as a function of 77 and Zvtx and the number of events is 
corrected as a function of reconstructed track multiplicity and Zvtx- The number of events without trigger 
or without reconstructed vertex is estimated from the simulation and included in the corrected number of 
events. Einally, the quantity dAch/d??, averaged over all events, is obtained for each 77 bin. The range of 
Zvtx contributing to the multiplicity varies with 77 (Eig.[T]l. Eor instance, at 77 = 2, tracks originate mostly 
from vertices in the range: -30 cm < Zvtx < -5 cm. Therefore, for each 77 bin, a Zvtx acceptance correction 
is applied. See ll55l for details of the procedure. 


6.2 Strangeness correction 


Since ALICE’S definition of primary charged particles excludes particles originating from the weak 
decays of strange particles, data have to be corrected for cases when daughter particles of such decays 
pass the track selection. Current Monte Carlo event generators have a strangeness content which differs 
from dafa by a factor approaching 2. Therefore, fhe sfrangeness confenf in fhe Monfe Carlo simulation 
was normalized fo dafa using AEICE’s K° and A measuremenfs in I 77 I < 0.9 [4], fhaf were exfrapolafed fo 
I 77 I < 2 using fhe shape from simulafion. The rafios of sfrangeness confenfs befween dafa and Monfe Carlo 
generators are slighfly cenfre-of-mass energy dependenf. Eor varying from 0.9 fo 8 TeV fhey increase 
from 1.6 to 1.85 according to PYTHIA 6 , and from 1.4 to 1.6 according to PHOJET. The uncerfainfy on 
fhese rafios coming from fhe uncerfainfy in AEICE measuremenfs of sfrange particle producfion [4], is 
esfimafed fo be 5%. The sfrangeness confaminafion is slighfly 77 dependenf, and varies from 1.7% af 77 = 0 
to 2.5% af 77 = 2 afV^ = 7 TeV. The sfrangeness correction is abouf 1%, has no significanf 77 variation in 
I 77 I < 2 and no significanf energy dependence befween ^/s = 0.9 and 8 TeV. This correction is explained in 
more defail in fhe Secfion 8.L3[ where fhe corresponding sysfemafic uncerfainfy is discussed. 


6.3 Event class normalization 

The final correcfion applied fo fhe dafa is fhe normalizafion fo one of fhe fhree evenf classes defined in 
fhis sfudy: NSD, INEE and INEE>0. In fhe normalization to NSD, correcfions have fo be made for fhe 
fraclion of SD evenfs remaining in fhe selection and for fhe fracfion of double-diffraction (DD) evenfs nof 
included in fhe selecfion. In fhe normalizafion of resulfs fo fhe INEE evenf class, corrections have fo be 
made for fhe fracfion of single- and double-diffracfive evenfs nof included in fhe selecfion. The INEE>0 
class is of inferesf because if minimizes diffracfive correcfions. In addifion, AEICE measuremenfs of SD 
and DD cross-sections BTlI reduced fhe sysfemafic uncerfainfies coming from diffraction. Correcfions for 
higher order diffracfive processes associafed wifh evenfs wifh fwo or more pseudorapidify gaps (regions 
devoid of parficles) are neglecfed in fhe normalizafion fo INEE, NSD and INEE>0 classes, as fheir 
confribufion fo inelasfic collisions is expecfed fo be smaller fhan 1 % lfT4ll5^ . Eurfhermore, such evenfs 
fend fo have a high frigger efficiency, which makes corresponding correcfions even smaller. 

yfs (TeV) MB and MBqr 

0.9 0.94+0'°2 o.91+o°3 

2.76 0.93«“ 0.88;|;« 

7 o.93;||:” o.85;S“ 

8 0.93«;“ 0.85^;“ 

Table 2: MB and trigger efficiencies for NSD events and MBqr trigger efficiencies for inelastic events at four 
centre-of-mass energies, obtained from diffraction-tuned versions of PYTHIA6 PemgiaO ll48ll and PHOJET lISTl . 
Uncertainties listed are total uncertainties. Statistical errors are negligible. The asymmetry of the MBqr errors is 
due to the asymmetric uncertainties in the diffraction efficiencies. 
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To normalize measurements to a given event class, trigger biases must be corrected for, both at event and 
track levels. For the INEL and INEE>0 classes, the correction is straightforward using the MBqr trigger 
efficiency (Table |^. 

For the NSD event class, contamination of the event sample by SD events must be taken into account. 
The measured quantity may be re-written as: 


1 


d(liMrk)MBi 


(■^ev)MB4 


drj 


1 

[d 

(2«,r) 

MB and , 

(sNf) 
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MB and 


MB and 


MB and 


dr] 


dr] 



(3) 


/■Class \ 


of a given class 


where is the number of tracks aggregated over all events (A' i 

V trjc /Trigger V /Trigger 

(superscript) selected with a given trigger type (subscript outside the parentheses). Given that 

MB and V /MBand^ 

tively, for SD or NSD events HTl . one obtains: 


MB and 

where s and cr are efficiencies and cross-sections, respec- 


1 

. ^MBand^^^ 1 1 d(ZMrk)MB,t.D ^MBand^^^ ^ ^trk )mb,nd ... 

I in ^ 

The coefficient in front of the single diffraction term in Eq. Q, varies from 0.04 dl^fs = 0.9 TeV to 
0.003 atV^ = 8 TeV. As the single diffraction term is not measured, but corresponds to a relatively small 
correction, this term was calculated using the simulation. The corresponding uncertainty was estimated 
by varying the single diffraction term conservatively between extreme cases, assuming either no SD, or 
assuming that all events are from SD. The last step consists of correcting for the MB and trigger efficiency 

to obtain the desired quantity, 1 /(a^^°) Idn. 

The DD event content of the MBqr and MB and data samples, is small, of the order of 5.5% and 4.5%, 
respectively. These fractions do not vary significantly between 0.9 and 8 TeV. The corrections for DD 
efficiency are included in the general efficiency correction. For the INEE and INEE>0 event classes, 
the MBqr trigger efficiency for DD events as a function of multiplicity is the same as for the other 
inelastic events. The MB and selection, which is used for the NSD event sample, has an efficiency for 
DD events that is lower than that of the other inelastic events. However, we checked in the simulation 
that the average efficiency correction for the NSD event class gives the same result as separate efficiency 
corrections implemented for DD and ND events. 


7 Multiplicity distributions of primary charged particles: analysis 

7.1 Unfolding multiplicity distributions 

The data samples used in these measurements are described in Table [T] The next step in the analysis 
consists of correcting the raw distributions for detector acceptance and efficiencies, using an unfolding 
method. 

The unfolding procedure follows the same approach as in [SJ, i.e. the corrected distribution is constructed 
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by finding the vector U, which minimizes given by 


2 

X = 


-Z 


+fixF{U) 


(5) 


where M represents the raw multiplicity distribution vector with uncertainty vector 5 , U the unfolded 
multiplicity distribution vector, and R the detector response matrix. Indices m and t run from 0 to the 
maximum number of multiplicity bins, in raw and corrected distributions respectively. The regularization 
term ySxF(f/) is used to decrease the sensitivity of the unfolding to statistical fluctuations. For F{U) a 
usual Tikhonov-type of function Il56l was used, which has a smoothing effect on the unfolded distribution 


N-l 




t=i 


Ut 


( 6 ) 


where N is the number of unfolded multiplicity bins, evaluated with the help of the response matrix, 
from the maximum raw multiplicity. 

The weight /3 (Table was chosen to minimize the mean squared error 1561. The solution is found to 
be stable over a broad range of /3 values (±50%), and the correct minimum was ensured in each case 
by scanning p over few orders of magnitude. The particular values of optimal weights depend on many 
features of the unfolding problem, such as distribution size, a pattern of fluctuations in the input raw 
data, properties of the response matrix and the regularization term. The most obvious dependence was 
eliminated by factorizing N in Eq. Q. 


(TeV) 


A;, 


0.5 

1 

1.5 

0.9 

10 

10 

10^ 

2.76 

10 

102 

102 

7 

10 

10^ 

10^ 

8 

50 

10^ 

10^ 


Table 3: Values of the weight parameter p used in the regularization term (Eq. Q), for each centre-of-mass energy 
and for each pseudorapidity range. 


For each generated multiplicity bin A®®" = t, the response matrix column Rmt consists of the distribution 
of the probability to measure multiplicity Ach = m. To extend the response matrix to the highest multi¬ 
plicities encountered in this study, beyond the reach of the available simulation, probability distributions 
were parameterized and extrapolated towards high multiplicities (Fig.|^. In the low-A®*^" region (A®*^" < 
10 to 20, depending on the rj range) the response matrix was taken directly from the simulation. In the 
large A®“ region (A®®" > 10 to 20), the column R^t is well described by a Gaussian distribution and 
mean values follow a linear trend (Fig. [^. Widths were parameterized using two different functions, a 
Fade function and a power law 


l+Cjl 

W{t) = Co + CiF 


Fade 

Fower law 


(V) 

( 8 ) 


Co, Cl, C 2 , C 3 and y are constants to be fitted. These functions have different asymptotic behaviours 
(Fig.|^, however, using either function makes a difference only for multiplicifies above 100 (in \r]\ < 1.5). 


The switch to parameterization occurs at A®“ = 10, 15 and 20, for \t]\ < 0.5, 1 and 1.5, respectively, for 
all energies. These values ensure that using the parameterized response matrix introduces no distortions 
in the low multiplicity region. 
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Generated ([ri| < 1.5) Generated (|ti| < 1.5) 


Fig. 5: Example of Gaussian parameterization of the response matrix, at Vi = 7 TeV, for It/I < 1; (left) paramete¬ 
rization of the mean values, with a linear function (red dashed line); (right) parameterization of the widths, with a 
Fade function (red solid line) and a power law function (blue dashed line). The bottom parts of the hgures show 
the ratios between data and hts. 


The range of multiplicities in the final unfolded distribution was further restricted by requiring that the 
bias (an estimate of how far is the result from the true solution ll56l ) is less than 10% in each bin. As 
unfolding is performed for each correction scenario (see Section [O] on systematic uncertainties), in the 
end the multiplicity range is limited by the unfolding resulting in the shortest range. The quality of the 
unfolding was verified by comparing fhe raw disfribufion M fo fhe producf /? ® (7. 

7.2 Event class normalization 

After the unfolding step, distributions have to be corrected for event selection efficiency (including trigger 
efficiency and vertex reconstruction efficiency). 

For the INEL and INEE>0 event classes this is straightforward, given that 

Mm = 'y^StRmtUt (9) 

t 

where St is the selection efficiency for true multiplicity t. Thus the unfolded distribution can be normal¬ 
ized to the inelastic event class by dividing the contents of each multiplicity bin by the corresponding 
efficiency. 

For the NSD event class the procedure is different as the unfolded distribution, U*, includes a contami¬ 
nation by SD events 

L; = ef (10) 

where upper indices denote the event class. The overall fraction of SD in inelastic collisions (a?) was 
measured by the AEICE Collaboration ll^ 

The desired unfolded distribution normalized to the NSD event class is obtained by combining Eqs 
and ( [TT| ) 

U* 

ttNSD _ _ 

.J^SD I a, „SD 
\-a,t 


.([Tg 

( 12 ) 
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8 Study of systematic uncertainties 

8.1 Common sources of systematic uncertainties 

8.1.1 Material budget 

The material budget in the ALICE central barrel was checked in the range \r]\ < 0.9, by comparing mea¬ 
sured and simulated gamma conversion maps. The conclusion is that, in this pseudorapidity range, the 
material budget is known with a precision of 5% Il46]l . The corresponding systematic uncertainty was 
obtained by varying the material budget in the simulation, conservatively over the whole pseudorapidity 
range by +10%, which induces a variation of dAch/d77 of ±0.2% for all event classes. Eor multiplicity 
distributions, the rj range considered does not exceed ±1.5, making the effect of the higher material bud¬ 
get uncertainty outside \r]\ < 0.9 relatively small. The systematic uncertainty from material budget is 
negligible compared to other sources of uncertainty. 

8.1.2 Magnetic field 

The magnetic field map was measured and simulated with finite precision. To check the sensitivity of 
the detector response to the precision of the simulation of the ALICE solenoidal magnetic field, data 
samples collected at = 7 TeV with opposite polarities were compared. The differences in dAch/ d?/ 
values are consistent with observed fluctuations between runs within the data-taking period at this energy. 
Therefore, the contribution from systematic uncertainties associated with the magnetic field are smaller 
than, or of the order of, the run-to-run fluctuations, and have been neglected. 

8.1.3 Strangeness correction and particle composition 

The main sources of uncertainty associated to the correction for strange particles originate from: (a) the 
difference (< 5%) in K® and A detection efficiency between data and simulation [4|; (b) the difference 
in pt distributions of strange particles in data compared to simulation, which implies a difference in the 
fractions of daughter particles meeting the vertex association condition; and (c) the uncertainty in the 
simulation of the strange particle content. 

Eor dNchl d ?7 measurements, the systematic uncertainty from strangeness correction is found to have a 
small T] variation, it is slightly larger at rj = 0 compared to \r]\ = 2. These uncertainties have a small 
energy dependence, they increase slightly with increasing V^, from 0.14% at = 0.9 TeV to 0.16% at 
^/s = S TeV. The uncertainties at 77 = 0 are listed in Table 

Eor multiplicity distributions, the strangeness contamination was studied with the Monte-Carlo simu¬ 
lation by evaluating the survival probability of strange particle decay products for the track selection 
used in the analysis. The probability that a track from strange particle decay passes the track require¬ 
ments is less than 0.1 % on average, leading to a negligible contribution to the uncertainty on multiplicity 
distributions. 

The particle composition affects the efficiency estimate, because different particle species have different 
efficiencies and effective pj cut off. The influence of the uncertainty in particle composition was esti¬ 
mated by varying, in the simulation, the relative fractions of charged kaons and protons with respect to 
charged pions by ±30%, which covers conservatively the uncertainties in the measured particle composi¬ 
tion at the LHC |[57l . and found to range, for dNchldij, between 0.1% for the INEL events class and 0.2% 
for the NSD and INEL>0 event classes (Table |^. The effect is negligible for multiplicity distributions. 

8.1.4 Detector Simulation 

The systematic uncertainty related to the limited precision with which the detector performance is sim¬ 
ulated was evaluated by varying the threshold on parameters used to select the various types of tracks. 
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Event 


^/s (TeV) 




class 

0.9 

2.76 


8 


7 

Material budget 

All 

0.2 

0.2 

0.2 

0.2 

Strangeness corrections 

All 

0.14 

0.17 

0.16 

0.16 

Pj uncertainties 

All 

-1-1.0 

-Hl.O 

-Hl.O 

H-l.O 

-0.5 

-0.5 

-0.5 

-0.5 


INEL 

0.1 

0.1 

0.1 

0.1 

Particle composition 

NSD 

0.2 

0.2 

0.2 

0.2 


INEL>0 

0.2 

0.2 

0.2 

0.2 


INEL 

h-3.5 

-h6.7 

-h7.3 

h-7.3 


-1.2 

-4.0 

-3.5 

-3.5 

Event generator dependence 

NSD 

h-4.0 

-1-3.9 

-h2.2 

h-2.2 


-2.0 

-3.0 

-2.0 

-2.0 


INEL>0 

h-0.5 

-hO.5 

-hO.6 

-1-0.6 


-0.5 

-0.5 

-0.6 

-0.6 

SPD simulation 

All 

0.8 

1.1 

0.6 

0.7 


INEL 

h-3.8 

-1-6.9 

-h7.4 

-3.6 

h-7.4 

-3.6 


-1.6 

-4.2 

Total systematic uncertainty 

NSD 

h-4.2 

h-4.2 

-h2.6 

h-2.6 


-2.3 

-3.3 

-2.2 

-2.3 


INEL>0 

-1-1.5 

+ 1.7 

-hI.5 

-1-1.5 


-1.2 

-1.4 

-1.2 

-1.2 


INEL 

0.4 

0.2 

0.4 

0.4 

Run-to-run fluctuations 

NSD 

0.4 

0.3 

0.3 

0.3 


INEL>0 

0.3 

0.2 

0.4 

0.4 


Table 4: Relative contributions, in percent, to systematic uncertainties in the measurement of dNch/drj, atri-0, 
for centre-of-mass energies considered in this study and for the three event classes defined in the text; the pj un¬ 
certainties combine the contributions from the difference between MC models and data with the uncertainty on 
the pt distribution below 50 MeV/c. Run-to-run fluctuation contributions indicated here for comparison are not 
included in total uncertainties. 

over a range obtained from the observed difference in the distributions of these parameters between sim¬ 
ulation and data. For dAch/d ?7 measurements based on tracklets, the;^f^ cut was varied between 1.3 and 
4. The spread of dAch/d?/ values over the range of z-positions of the vertex covered by a given rj bin was 
used as a measure of the bias introduced by the z-dependence of the tracklet reconstruction efficiency. 
The corresponding uncertainty ranges from 0.6% to 1.1% depending on the data sample (Table|^. 

Since for multiplicity distribution measurements all three track-counting algorithms are used, track pa¬ 
rameters were also varied for global tracks, and for ITS-only tracks (details may be found in 1461 1: 

- for global tracks, the minimum number of TPC clusters was varied between 60 and 80 and the;^f^ cut 
between 3 and 5. 

- for ITS-only tracks, the cut was varied between 2 and 3. 

The uncertainty on the DCA distribution is not considered here, as it is included implicitly in the uncer¬ 
tainty on the strangeness content, which affects the DCA distribution. 

We find fhaf, in fhe case of mulfiplicify disfribufions, uncerfainfies in fhe defector simulafion are negligi- 
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ble compared to other sources of uncertainties. 

8.1.5 Model dependence 

The remaining SD fraction in the sample selected with the MB and trigger in view of the normalization 
to the NSD event class is 3%, 1% and negligible = 0.9, 2.76 and > 7 TeV, respectively. Uncertain¬ 
ties coming from diffraction contributions are included in the trigger efficiency uncerfainfies (Table 
obfained in BTl . for < 7 TeV. = % TeV, fhe efficiency values were faken fo be fhe same as for 
yfs -1 TeV. In addifion, fhe model uncertainties at^fs = 0.9, 2.76, 7 and 8 TeV were obtained from the 
difference befween PYTHIA6 PerugiaO and PHOJET. A tesf of fhe efficiency evaluation was obtained 
by comparing simulated MB and to MBqr trigger efficiency ratios to the measured values. Excellent 
agreement was found at all energies WJ\ . 

Eor multiplicity distributions, the systematic uncertainty from the model dependence is included in both 
efficiency correction and px dependence uncertainties, as different event generators and tunes are used 
to estimate independently efficiencies and response mafrices (see Section [8^ . 

8.1.6 p'x dependence 

None of fhe MC generafors used in fhe defector simulafion reproduces correctly the p^ distribution of 
charged particles observed in the data If5^l59l . This introduces an uncertainty in the determination of the 
detector response, as it is integrated over transverse momenta, and the probability of detecting a particle 
decreases with decreasing pj- This affects in particular, together with the uncertainty on the material 
budget and the magnetic field, extrapolations of measurements to px = 0. 

In order to study px spectrum effects, two different tunes of the PYTHIA6 event generator were used, 
ATEAS CSC and PerugiaO, which give an average px versus charged-particle multiplicity respectively 
below and above the data (Pig. (left)), for most of the multiplicity range (Vch > 2). The difference 
befween measurements obtained with response matrices corresponding to each of these Monte Carlo 
generators, is used as the corresponding systematic error contribution. Pigure (right) shows that this 
procedure introduces an uncertainty which is weakly dependent on p, and amounts to ±0.3%, when 
averaged over |p| < 2. 

Eor undetected particles, below a threshold of about 50 MeV/c, a value chosen to coincide with a track 
detection efficiency of 50%, fhe corresponding sysfematic uncerfainfy is obtained by varying their frac¬ 
tion by a conservative amount (-50% and ±100%). (The fraction of the px spectrum below 50 MeV/c is 
about 1 % of the total, for both PYTHIA6 and PHOJET.) The resulting systematic uncertainty on dVch/ dp 
is p dependent, and found to range from -0.5% and ±1.0%, at p = 0, to -0.75% and ±1.5%, at |p| = 2. 

The systematic uncertainty induced on P(Vch) by the difference in px befween dafa and simulafion is 
slightly sensitive to the tune of PYTHIA6 considered: for instance, at Vch = 90, varying the px spectrum 
below 50 MeV/c by -50% and ±100% induces a change of -5% and ±9%, respectively for the ATEAS 
CSC tune and of -4% and ±8%, respectively, for the PerugiaO tune. 

8.2 Systematic uncertainties in unfolding of multiplicity distributions 
8.2.1 Uncertainty evaluation 

The results of the unfolding procedure and unfolding uncertainty estimate were cross-checked, in stan¬ 
dard ways: 

- changing the regularization term by either varyingyS or the function F(Uy, 

- using two alternative unfolding procedures: Bayesian Il60l and singular value decomposition llFD . 
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Fig. 6: (left) At = 7 TeV, average raw pj vs. raw charged track multiplicity in |? 7 l < 1, for data (black circles), 
PYTHIA 6 PerugiaO (blue dashed line) and ATLAS-CSC (red dotted line). The bottom part of the figure shows the 
ratios of the two simulated distributions to the data; (right) Comparison of dAch/ d ?7 evaluations, as a function of 77 , 
using correction maps obtained with (a) PYTHIA 6 PerugiaO (blue dashed line), (b) PYTHIA 6 ATLAS-CSC (red 
dotted line) and (c) the average between (a) and (b); The bottom part of the figure shows the ratios of corrected 
distributions from (a) and (b) to the average. 

The changes in unfolded distributions due to variations of the different elements used in the unfolding 
procedure, each having their own systematic uncertainty, were studied by considering all possible combi¬ 
nations: counting algorithm (3 cases), event generator used for efficiency correction (2 cases: PYTHIA 6 
and PHOJET funed for diffraction), evenf generator used for the response matrix (2 cases PYTHIA 6 CSC 
and PYTHIA 6 PerugiaO), low pj spectrum extrapolation below 50 MeV/c (3 cases: varying the integral 
below 50 MeV/c by -50%, 0% and -i-100%), response matrix parameterization (2 cases: Power law, and 
Pade), which correspond to 72 separate measurements that are correlated as a result of the unfolding 
procedure. To take these correlations into account, for a given energy and a given rj range, the systematic 
uncertainty from all the sources considered was estimated as the overall spread between the resulting 
distributions. The average between these distributions (center of the band covered by the 72 curves) was 
used as the measurement. The total uncertainty originating from the unfolding procedure (evaluated for 
each multiplicity bin as a linear sum of the unfolding bias |[5^ and the unfolded distribution covariance 
calculated from the statistical uncertainty of the raw distribution) was added linearly to the systematic 
uncertainty defined as half fhe size of fhe band. This fakes info account the fact that the unfolded dis¬ 
tributions come from the same raw multiplicity histogram, hence the raw data statistical fluctuations are 
propagated to each of the unfolded distributions in a similar way and affect the spread of distributions 
uniformly. The resulting systematic uncertainty is expected to be highly correlated with multiplicity. 

The systematic uncertainties from material budget, tracklet and track selection, detector alignment (eval¬ 
uated by changing the geometry within alignment uncertainties), particle composition, strangeness cor¬ 
rections are found to be negligible. Among the non-negligible contributions to multiplicity distribution 
systematic uncertainties, the contribution from selection efficiency (Secfion [8.2.1[ ) can be evaluated sep¬ 
arately as it is a multiplicative correction. All the other contributions, from the uncertainty on {pj), the 
extrapolation down to pj = 0 , the counting algorithms, and the model dependence including the con¬ 
tribution from diffraction, are all mixed together through the procedure described above. Eor NSD and 
INEE event classes, the diffraction contribution is mainly significant in the zero multiplicity bin, which 
is absent for the INEE>0 event class. 
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Source of 

Event 


yfs (TeV) 


uncertainty 

class 

0.9 

2.76 

7 

8 

S' 

H<0.5 

INEL 

NSD 

INEL>0 

0.96-0.7-0.5 

4.68-3.5-2.4 

--0.7-0.5 

3.05-2.1 - 1.2 

9.20-6.3-3.7 
--2.1 - 1.2 

2.11 - 1.5-0.9 
4.55-3.2-1.9 

-- 1.5-0.9 

2.52-1.7-1.0 
13.45-9.1-5.3 
-- 1.7- 1.0 

Efflcienc 

icertainty 


INEL 

NSD 

INEL>0 

0.14-0.2-0.2 
0.42-0.8-1.2 
--0.2-0.2 

0.49 - 0.5 - 0.5 
1.66-2.1-2.3 
--0.5-0.5 

0.18-0.2-0.3 
0.31-0.8-1.1 
--0.2-0.3 

0.02-0.2-0.3 
2.91-3.1-3.1 
--0.2-0.3 

D 

H < 1.5 

INEL 

NSD 

INEL>0 

0.02-0.2-0.3 
0.10-0.4-0.7 
--0.2-0.3 

0.03-0.1-0.2 
0.12-1.2-2.0 
--0.1 -0.2 

0.03-0.3-0.4 

0.01-0.4-0.7 

--0.3-0.4 

0.03 - 0.2 - 0.4 
0.19-1.6-2.8 
--0.2-0.4 


Table 5: Systematic uncertainties, in percent, due to efficiency corrections including trigger and reconstmction 
efficiency, event generator dependence and diffraction, for INEL, NSD and INEL>0 event classes, in pseudora¬ 
pidity intervals \ri\ < 0.5, 1.0 and 1.5, - 0.9, 2.76, 7 and 8 TeV. Numbers are given, from left to right, for 

multiplicities 0, 1 and 2. Eor the INEL>0 event class, the efficiencies are the same as for the INEL class, except 
that there is no bin with multiplicity 0. Note that the uncertainty listed here for bin 0 is not the total uncertainty in 
that bin. Bin 0 values are recalculated separately from simulation, which adds a significant spread of values that is 
seen in the total uncertainty estimate. 


Source of 

Event 


s/s (TeV) 


uncertainty 

class 

0.9 

2.76 

7 

8 



INEL 

4-4-16 

2-2-6 

1-3-4 

22-9-13 


|77l<0.5 

NSD 

16-6-28 

28-9-15 

19-7-9 

35-8-12 

E 

t? .B 

t/} -f- 


INEL>0 

8-3-14 

3-2-6 

3-4-4 

3-1-4 


INEL 

7-7-21 

3-3-8 

1-3-5 

29-11-13 

\ri\ < 1 

NSD 

38-10-32 

52-8-15 

38-8-11 

59-8-14 

« *- 


INEL>0 

14-6-20 

7-3-8 

10-3-5 

4-2-4 








H =s 


INEL 

11-9-27 

3-7-12 

3-5-8 

34-9-14 


l7l<l-5 

NSD 

62-8-34 

76-12-24 

57-8-14 

80-9-17 



INEL>0 

12-7-27 

7-7-11 

10-6-8 

8-1-6 


Table 6: Total systematic uncertainties in charged-particle multiplicity distribution measurements, in percent, for 
INEL, NSD and INEL>0 event classes, in pseudorapidity intervals \t]\ < 0.5, 1.0 and 1.5, at Vs = 0.9, 2.76, 7 and 
8 TeV. Values for INEL and NSD event classes are given for three characteristic multiplicities, from left to right: 
zero multiplicity - mean multiplicity - five times the mean multiplicity. Eor the INEL>0 class, the first value is for 
Nch = 1, as there is no entry for N^h - 0. 

8.2.2 Bin-to-bin correlations in systematic uncertainty 

Due to specific nature of the correction procedure, the final multiplicity distributions contain bin-to-bin 
correlations coming from various sources with different properties. These sources can be categorized by 
their effect: 


- statistical correlations, resulting from the propagation of the raw distribution statistical uncertainties 
through unfolding process, their characteristics largely depend on the response matrix structure; 

- fully correlated shift of the distribution as a result of the uncertainty in normalization; 

- long-range correlations in systematic uncertainties as a result of multiplicity scale change (that is 
determined by the position of the response matrix bulk relative to the diagonal, see Section [TT] ) in the 
unfolded distribution. 
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We found that first category correlations (statistical) are negligible compared to the last category (scal¬ 
ing), while the second category correlations can be easily factorized for fitting the multiplicity distribu¬ 
tion. The three main sources of correlated uncertainty, falling into the third category, are the change of 
counting algorithm, the change of event generator tune used to produce the response matrix (see Fig. 
and the variation of pT distribution under the detection threshold. In order to evaluate the effect on the 
final distributions we construct 18 intermediate distributions corresponding to all possible combinations 
of counting algorithm, event generator tunes and variations of px spectra. All other sources of systemat- 
ics are included (except the uncertainty corresponding to bin Nch = 0 renormalization, as it is not applied) 
with the same procedure described for the final distribution previously in this section. These 18 distri¬ 
butions are then treated independently in Section [93] to evaluate the effect of correlated uncertainties on 
NBD fits. 


8.3 Summary of systematic uncertainties 
8.3.1 Pseudorapidity density 

The various contributions to systematic uncertainties in dAch/dp are summarized in Tablej^for the three 
event classes and the four centre-of-mass energies studied in this publication. For the INEF>0 event 
class, a precision of 1.5% is achieved, as the sensitivity to diffraction is negligible and the p range is 
reduced in the definition of this event class. 

In the p range covered in this study (|p| < 2), we find that systematic uncertainties show essentially no p 
variation (Fig.|^, and are therefore strongly correlated bin-to-bin. 


8.3.2 Multiplicity distributions 

The efficiency uncertainties (Section [8.2.1[ ) are only relevant at low multiplicity. For multiplicities above 
8 to 9, the efficiency reaches 100% and the corresponding systematic uncertainty becomes negligible. 
Therefore, efficiency uncertainties are only given for a few characteristic low multiplicities. The total 
systematic uncertainties vary with multiplicity, therefore they are given in Section 8.2.1| only for a few 
characteristic values of the multiplicity. 


8.3.3 Consistency checks 

In the measurement of dAch/dp, statistical errors are negligible. Therefore, the study of run-to-run 
fluctuations (measured RMS of dAch/dp results in different runs) provides a check that all run dependent 
corrections are properly handled. 

For the INEF and NSD event classes, contributions from run-to-run fluctuations are significantly smaller 
than the total systematic error (Fig.[^. For the INEF>0 event class, for which the precision is highest, the 
relative importance of run-to-run fluctuations is larger than for the INEF and NSD event classes (Fig.|^, 
but reaches at most 5% of the total systematic uncertainty. 

As data correction procedures are significantly different between dAch/ dp and multiplicity distribution 
measurements, a test was performed to verify the consistency of the two measurements. At the four 
centre-of-mass energies and for the three pseudorapidity intervals used in this study, integrals of the 
multiplicity distributions were found to be consistent with the direct measurements of dAch/dp, within 
errors. 
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Fig. 7: Total relative systematic uncertainty on dAch/d?/ (thick black lines), as a function of pseudorapidity, com¬ 
pared to run-to-run fluctuations (thin red dashed lines) at Vi = 0.9 TeV (top row), 2.76 TeV (second row), 7 TeV 
(third row) and 8 TeV (bottom row), for the INEL, NSD and INEL>0 event classes, as indicated. 


9 Experimental results 

9.1 Pseudorapidity density of primary charged particles: measurements 

In their common rj range, I 77 I < 0.9, the three track counting algorithms discussed in Section |5^ achieve a 
similar precision of 1%, and were found to give consistent results. The main difference is that for ITS-i- 
and ITSTPC- 1 -, there is a detector calibration contribution to systematics for the TPC and the SDD, not 
present for the SPD. To achieve the largest possible 77 range, in the measurement of dAch/d?/ versus ij, 
the Tracklet algorithm is used alone. 

The measurement at V^ = 0.9 TeV is shown in Fig. [^compared with previous results. At \t]\ > 0.9 the 
measurement for the INEL event class is slightly lower than in ALICE’S previous publication Q- The 
difference comes mainly from: (a) the tuning of the MC generators for diffraction, as larger pseudo¬ 
rapidities are more sensitive to SD; (b) the subtraction of particles coming from the decay of strange 
particles was improved, using ALICE’S measurement of strangeness ||4l; and (c) the improvement of the 
T] dependence of the Tracklet algorithm. 

The discrepancy with UA5 for the INEL event class at large rj could perhaps be related to the fact that 
UA5 used a 1/Mx variation of the single-diffractive cross section (see ll53l ). Note also that UA5 data 
seem to be internally inconsistent (see discussion in Il62ll ). The measurement at V^ = 2.76 TeV is shown 
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Fig. 8: dNch/dr] vs. rj at^Js - 0.9 TeV, for the three normalizations defined in the text, and a comparison with 
ALICE previous measurements |l2]|3], UA5 ll^ and CMS ll^ . Note that to avoid overlap of data points on the 
figure, the INEL>0 data were displaced vertically, and for these data the scale is to be read off the right-hand side 


vertical axis. Systematic uncertainties on previous data are shown as error bars (except for UA5, with coloured 
bands), while they are shown as grey bands for the data from this publication. 


in Fig.j^top, and found to be consistent with the ALICE measurement at = 2.36 TeV |j2jl, as expected 
because of the small change in center-of-mass energy. The new measurements of dAch/d??, at = 7 TeV 
(Eig. [^middle), show agreement both with the previous AEICE results for INEE>0 lO, and with CMS 
NSD data |[33ll . The measurements of dAch/d?? at-^ = 8 TeV (Eig. [^bottom) show the 3% increase with 
respect to the 7 TeV data, which corresponds to what is obtained in the extrapolation from lower energy 
data. Comparisons of the rj distributions at the four centre-of-mass energies (Eig. [TOj l, for the three events 
classes, show no significant change of shape and a smooth increase of the charged-particle density with 
increasing energy. 


The data for the INEE event class at^[s = 0.9 and 7 TeV were compared to simulations with current event 
generators (Eig. 111. At= 0.9 TeV, EPOS EHC |[63l and PYTHIAS 4C Il64ll65ll are consistent with 
the data. PHOJET overestimates the data, while PYTHIA6 PerugiaO and Perugia 2011 underestimate 
the data. At = 7 TeV, EPOS EHC, PHOJET and PYTHIA6 Perugia 2011 are consistent with the 
data. PYTHIAS 4C overestimates the data, while PYTHIA6 PerugiaO underestimates the data. Note that 
PYTHIA6 Perugia 2011, PYTHIAS 4C and EPOS EHC were tuned using EHC data. 


9.2 Energy dependence of dAch/ d;/ at;/ = 0 


The traditional definition for dAch/d?; at ?/ = 0 is an integral of the data over the pseudorapidity range 
|;7l<0.5 

+ 0.5 


dAch 


d?/ 


,/=0 


- 0.5 


dAch 


d/7 


(13) 


The results of the measurements of dAch/d/y at /y = 0 are given in Table The energy dependence of 
dAch/ d/y at ?y = 0 is of interest not only because it provides information about the basic properties of pp 
collisions, but also because it is related to the average energy density achieved in the interaction of pro- 
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Fig. 9: dAch/d ?7 vs. rj measurements: ^fs - 2.76 TeV compared with-\/i = 2.36 TeV taken from ALICE ||2l (top); 


\[s - 1 TeV and comparison with CMS and ALICE 13 data (middle); V* = 8 TeV (bottom). Systematic 


uncertainties are shown as error bars for the previous data and as grey bands for the data from this publication. The 


scale is to be read off the right-hand side axis for INEL>0. 
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Fig. 10: Comparison of dNch/dj] vs. rj measurements between the various centre-of-mass energies considered in 


this study: INEL (left), NSD (middle), and INEL>0 (right). The lower parts of the figures show the ratios of data 
at energies indicated to the data at 0.9 TeV, with corresponding colours. Systematic uncertainties are indicated as 
coloured bands. 


^/s (TeV) 

INEL 

NSD 

INEL>0 

0.9 

9 Q4+0.11 

3 

■’•^^-0.16 

0 7C+O.O6 
'^-0.05 

2.36^^ 

7) 77+0.25 
/_0.12 

4 49+017 

— 

2.76 

o 7C+O.26 
'^•^'^-0.16 

^■'J-'-0.i9 

4 70+0.08 

7 

-0.17 

c 74+0.15 
-0.15 

5 98+0.09 

8 

H.OO-o 17 

5 90+015 

'^■^■’-0.08 


^Data taken from o 


Table 7: Summary of ALICE measurements of dAch/d?; at 77 = 0 (integral of the data over |? 7 | < 0.5), for centre-of- 
mass energies and event classes considered in this study. The errors shown are systematic errors. Statistical errors 
are negligible. 

tons, and constitutes a reference for the comparison with heavy ion collisions. At mid-rapidity, dAch/dty 
can be parameterized as dAch/dty ~ s^. Combining the ALICE data with other data at the EHC and at 
lower energies, we obtain 5 = 0.102 ± 0.003, 0.114 ±0.003 and 0.1 14±0.0015[^ for the INEE, NSD and 
INEE>0 event classes, respectively, to be compared to d - 0.15 for central Pb-Pb collisions Il54l . This 
is clear evidence that the particle pseudorapidity density increases faster with energy in Pb-Pb collisions 
than in pp collisions. Pits are shown on Pig. and Table gives extrapolations to centre-of-mass ener¬ 
gies of 13 and 14 TeV (EHC design energy). While this paper was being prepared, the first measurement 
at 13 TeV by CMS appeared ICTl . resulting in dVch/d? 7 ||^|<o .5 - 5.49 ±0.01 (stat)±0.17 (syst) for inelastic 
events, which is consistent with our extrapolation of 5.30 ±0.24. Over the EHC energy range, from 0.9 
to 14 TeV, while the centre-of-mass energy increases by a factor 15.5, extrapolation of present data for 
dNchldrj at ?7 = 0 shows an increase by factors 1.75 ±0.03, 1.87 ±0.03 and 1.87 ±0.01, respectively for 
the three event classes. The multiplicity increase is similar for NSD and INEE>0 classes but slightly 
lower for the INEE class. 

9.3 Multiplicity distributions of primary charged particles: measurements 

The results of ALICE measurements of multiplicity distributions of charged primary particles are dis¬ 
played as probability distributions (P(Nch)) in Pigures [T3] (INEE), [T4| (NSD) and[T5](INEE>0). Por the 

®The uncertainty on d was obtained, assuming that the data errors are independent at different centre-of-mass energies, 
which is not strictly the case. 
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Fig. 11: Comparison with models of ALICE measurements of dNch/dt] versus rj for the INEL event class, at 
^[s - 0.9 (left) and 7 TeV (right): ALICE data (black circles with grey band), PYTHIA6 tune PerugiaO ||4^ 
(red continuous line), PHOJET 1511 (blue dot-dashed line), PYTHIA6 tune Perugia 2011 l48l (pink dashed line), 
PYTHIAS 4C 1641 |65l (green dashed line), EPOS LHC 1^ (long dashed light blue line). The lower parts of the 
hgures show ratios of data to simulation. Systematic uncertainties on ratios are indicated by coloured bands. 


Vs (TeV) 

INEL 

NSD 

INEL>0 

13 

5.30 ±0.24 

6.50 ±0.20 

6 . 86 ± 0.10 

13.5 

5.33 ±0.25 

6.56 ±0.20 

6.92 ±0.10 

14 

5.37 ±0.25 

6.62 ± 0.20 

6.98 ±0.10 


Table 8: Extrapolations of dVch/d;/, at t; = 0, for the three event classes, to higher energies at the LHC {sfs =13 
and 14 TeV), using the hts described in the text. 


first two event classes the measurements were obtained in three pseudorapidity intervals |? 7 | < 0.5, 1 and 
1.5, and for INEL>0 in |? 7 | < 1. At yfs — 1 TeV, P (Vch) varies over 6 to 7 orders of magnitude and the 
multiplicity range reaches up to 160 in |? 7 | <1.5 for both INEE and NSD event classes. In I 77 I < 0.5 and 
|? 7 l < 1, the observed multiplicity reaches 10 times the mean multiplicity. It is expected that the average 
energy density in proton collisions at the EHC, at = 14 TeV, is about 5 to 15 times smaller than en¬ 
ergy densities reached in gold ions at RHIC fTfil . Therefore, in proton-proton collisions of multiplicity 
exceeding 10 times the average multiplicity, energy densities should overlap with those of heavy ion col¬ 
lisions at RHIC, allowing to compare properties of systems with very different collision volumes (two to 
three orders of magnitude) but the same energy density. Euture runs of the EHC should allow extending 
much further the range of multiplicities probed so far. 


The high-multiplicity tail of the distribution increases as expected with increasing energy (Eig. 16 1 . This 
behaviour is studied quantitatively in Section [V^ on KNO scaling and < 7 -moment analysis. 

The measurements presented in this publication are consistent with previous AEICE data, for INEE |12| at 
- 0.9 TeV and INEE>0 0 at - 7 TeV, in the multiplicity range where they overlap (Eig. 17 1 . The 
wavy structure already observed by AEICE in IJJO, for multiplicities above Vch - 25, and previously by 
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Fig. 12: Charged-particle pseudorapidity density in the pseudorapidity region \ri\ < 0.5 (dAch/d;; at 77 = 0 calculated 
as the integral of the data over |? 7 | < 0.5) for INEL, NSD, and INEL>0 collisions, as a function of the centre-of- 
mass energy. Lines indicate fits with a power-law dependence on ^[s. Grey bands represent the one standard 
deviation range. Data points at the same energy have been shifted horizontally for visibility. The LHC nominal 
centre-of-mass energy is indicated by a vertical line. Data other than from ALICE used in this figure are taken 
from references ll25l [33] [66| 16814751 . 

UA5 If27l . is still hardly significant, and it is not present in the raw data. This feature was also observed in 
a study of CMS data llTTl . However, applying the same procedure to Monte Carlo data produces similar 
structures in the unfolded distribution, while no oscillation was present at particle level. We conclude 
that this structure is an unfolding procedure artifact. The period of the structure is related to response 
matrix width. 


9.4 Comparison of multiplicity distributions with other experiments and models 


CMS data are available for the NSD normalization only. Ai^fs - 0.9 and 7 TeV, in the three pseudo¬ 
rapidity intervals where they could be compared, the ALICE multiplicity measurements are generally 
in agreement with CMS data |[34l (Eig. 17 1 . However, at = 0.9 TeV, above Nch » 40 in |? 7 | < 1 and 
above Nch ~ 50 in |? 7 | <1.5, the AEICE multiplicity distributions tend to be higher than CMS data. The 
difference is nof significanf in each isolafed bin, however fhe general frend seems fo be differenf. As fhe 
bin-fo-bin correlafions in CMS’s resulfs are unknown, if is very difficulf fo obfain reliable unbiased quan- 
fifafive comparison. A possible source of fhe discrepancy is fhe differenf freafmenf of single-diffracfive 
evenfs in fhe fwo analyses making fhe NSD evenf sample definitions nof sfricfly compafible. More pre¬ 
cisely, fhere was no diffraction tuning in the simulations used by CMS and, moreover, AEICE’s criteria 
for events to be considered single-diffractive (see 1471) include a fixed cut on diffractive mass that differs 
from fhe value fhaf CMS used. This explanation is supporfed by fhe facf fhaf fhe difference is insignifi- 
canf af = 7 TeV, where single-diffracfive confribufion is negligible. The comparison wifh CMS dafa is 
furfher discussed in fhe Section [93] where NBD fifs fo fhe dafa are reported. 


Measured mulfiplicify disfribufions are compared fo models (Eig. [T8| ), for fwo energies, fhe lowesf one 
(V^ - 0.9 TeV), and fhe energy af which fhere is fhe largesf evenf sample {yfs - 1 TeV), for fhe INEE 
evenf class, in fhe pseudorapidify range I?)! < 1. 
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uncertainties combined; (a) data at Vi = 0.9 TeV (top left); (b) data at V* = 2.76 TeV (top right); (c) data at 
Vs = 7 TeV (bottom left); (d) data at Vs = 8 TeV (bottom right). Ratios of data to the fits are also shown, with 
shaded areas representing combined systematic and statistical uncertainties. 
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uncertainties combined; (a) data at Vi = 0.9 TeV (top left); (b) data at V* = 2.76 TeV (top right); (c) data at 
Vs = 7 TeV (bottom left); (d) data at Vs = 8 TeV (bottom right). Ratios of data to the fits are also shown, with 
shaded areas representing combined systematic and statistical uncertainties. 
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solid lines show the single and double NBD fits (see Section [93| ). Shaded areas represent statistical and systematic 
uncertainties combined; (a) data at V* = 0.9 TeV (top left); (b) data at V* = 2.76 TeV (top right); (c) data at 
= 7 TeV (bottom left); (d) data at = 8 TeV (bottom right). Ratios of data to the fits are also shown, with 
shaded areas representing combined systematic and statistical uncertainties. 

At = 0.9 TeV, PYTHIA6 PerugiaO fails to reproduce the data. PYTHIA6 Perugia 2011, PYTHIAS 4C 
and EPOS EHC underestimate the multiplicity distribution at high multiplicities (Vch ^ 25). PHOJET 
provides the best match to the data over the observed multiplicity range, up to Vch » 60. 

At - 7 TeV, both PYTHIA6 PerugiaO and PHOJET fail to reproduce the data. They severely underes- 
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8 TeV), for INEL and NSD event classes and for \ri\ < 0.5 (top row) and It/I <1.5 (bottom row). 


timate the high multiplicity part of the distribution. PYTHIA6 Perugia 2011, EPOS EHC and PYTHIAS 
4C give a reasonable fit of the low multiplicity region but underestimate the data above Ach » 60. 


9.5 Parameterization of multiplicity distributions with NBDs 

Single Negative Binomial Distributions (NBD) have been traditionally used to parameterize particle 
multiplicity distributions in hadron collisions 


Pnbd {n,{n),k) 


Tin + k) 

(n) 

n 

k 

r{k)r{n + l) 

(n) + k 


{n) + k 


(14) 


where {n) is the average multiplicity and the variance is given by 

D^ = {n^)-{n)^^{n)+^^ (15) 

The parameter k is related to the two-particle correlation function, in the pseudorapidity interval consid¬ 
ered fTSll . In the limit k ^ oo, the NBD becomes a Poisson distribution. 

In previous AEICE data, for the NSD event class, no strong deviation from a single NBD fit was observed 
atV^ = 0.9 and 2.36 TeV, for \r]\ < 1, while a hint of a substructure appears at \r]\ < 1.3 El. AtV^ = 7 TeV, 
for the INEE>0 event class, the single NBD fit slightly underestimated the data at low multiplicity 
(Vch < 5), and slightly overestimated the data at high multiplicities (Vch > 55) El- 

In the present data, for all event classes, already at^/s = 0.9 TeV, there is a hint that single NBD fits start 
diverging from the data at the higher multiplicity, for \r]\ < 0.5 and 1.0. More significant departure from 
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and the INEL>0 event class at V* = 7 TeV (top right); in both cases, ratios between ALICE new data and ALICE 
previous data are also shown. The total uncertainties are shown as error bars for the previous data and as a band 
for the present measurement. Eor the NSD event class, comparison with ALICE previous publication ||2] and with 
CMS data at V* = 0.9 TeV ll34l (bottom left), and comparison with CMS data at V* = 7 TeV ll34ll (bottom right); 
ratios of the NBD fits of ALICE data taken without errors to CMS data, for the various tj intervals (indicated on the 
figure), are also shown. Error bars represent the contributions of the CMS errors to the ratios, the bands represent 
the ALICE total uncertainty assigned to the ratio of 1. 
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Fig. 18: Comparison with models of measured multiplicity distributions for the INEL event class in the pseudora¬ 
pidity range \ri\ < 1.0: ALICE data (black circles with grey bands), PYTHIA6 Il48l tune PemgiaO (red continuous 
line), PHOJET ISB (blue dot-dashed line), PYTHIA6 Peragia 2011 (pink dashed line), PYTHIAS 4C ll64ll^ 
(green dashed line), and EPOS LHC 16^ (long dashed light blue line). The shaded areas represent total uncertain¬ 
ties: comparison at 0.9 (left) and at 7 TeV (right). The ratios between measured values and model calculations are 
shown in the lower parts of the figures with the same convention. 

single NBD starts at = 2.76 TeV. The values of x^j dof reflect this, as they are rather high and increase 
further with increasing centre-of-mass energies. 

The appearance of substructures in multiplicity distributions attributed to the occurrence of several 
sources in the process of particle production |[79] - [82l . can be parameterized by fitting the data with two 
NBDs. Indeed, a much better fit to the data is obtained by using a weighted sum of two NBD functions 

P{n) = T[aENBD(«,(«>! ,^i) + (l -a)ENBD(«,<«> 2 .^ 2 )] (16) 


This type of function, however, is not meant to describe the value P(0) for INEL and NSD distributions, 
which occurs when the rj acceptance is limited, therefore the bin n-0 was excluded from the fit and an 
overall normalization factor (T) was introduced, as a free parameter, to account for this. Kt^[s -1 TeV, 
P. Gosh et al. in fJl\ perform a similar fit for CMS data, however without using an overall scale factor 
(d) and with no account for the effect of bin-to-bin correlations. 


As was described in Section 8.2.2[ 18 multiplicity distributions, corresponding to main sources of cor¬ 
related uncertainty, were fitted independently. Eor each of these distributions systematic uncertainties 
corresponding to the non-leading sources were found to have negligible correlations, as expected. Thus 
these residual systematic uncertainties were added to the diagonals of the statistical covariance matrices 
for these distributions, that were produced in the unfolding procedure. In Tables [9| to [TT] we present the 
central values of the fit parameters, i.e. center of the maximum spread of each parameter between the 18 
independent fits. Half of the spread is used to estimate the systematic uncertainty of the parameters (the 
statistical uncertainty being negligible). However these estimates should be treated with caution due to 
the significant correlations between parameters. The multiplicity distributions for each of the 18 cases 
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differ by their slope and intersect at low-to-middle values of multiplicity {Nch ~ 10 to 30 depending on 77 
window and energy). There is no single point of intersection for a given set, but rather a small interval. 
Therefore extreme cases were selected as the curves with highest and lowest values at high multiplicities 
(and consequently lowest and highest values at low multiplicity). To indicate the change in shape of the 
parametrization allowed by the correlations present in systematic uncertainties, the parameters for these 
extreme cases are given in Tables [9| to [TT] as well. 


VS 

(TeV) 

A;, 

a 

<n>i 

ki 

<n>2 

k2 

4 



0.427 

1.87 

2.59 

4.92 

2.52 

0.805 


0.5 

1 0.469 ±0.118 

2.08 ±0.30 

3.09 ±0.67 

5.56±0.82 

2.97 ±0.60 

0.793 ±0.013 1 



0.587 

2.37 

2.64 

6.38 

3.48 

0.788 

0.9 


0.615 

4.46 

2.38 

11.98 

3.76 

0.813 

1 

1 0.627 ±0.048 

4.77 ±0.45 

2.45 ±0.31 

12.93 ±1.20 

3.95 ±0.38 

0.805 ±0.012 1 



0.612 

5.22 

2.49 

14.13 

4.33 

0.798 



0.444 

6.18 

3.54 

15.63 

3.04 

0.793 


1.5 

1 0.502 ±0.069 

6.81 ±0.66 

3.22 ±0.32 

17.36±1.74 

3.34 ±0.40 

0.794 ±0.008 1 



0.553 

7.25 

2.90 

18.79 

3.60 

0.791 



0.546 

2.39 

2.34 

7.46 

3.12 

0.799 


0.5 

1 0.588 ±0.052 

2.54±0.15 

2.20 ±0.36 

7.87 ±0.41 

3.32 ±0.24 

0.800 ±0.009 1 



0.578 

2.54 

2.13 

7.89 

3.24 

0.799 

2.76 


0.574 

5.05 

2.40 

15.54 

3.53 

0.802 

1 

1 0.591 ±0.025 

5.20 ±0.20 

2.26±0.18 

16.16±0.63 

3.63±0.19 

0.805 ±0.006 1 



0.591 

5.15 

2.08 

16.07 

3.52 

0.811 



0.580 

7.86 

2.55 

23.04 

3.69 

0.797 


1.5 

1 0.601 ±0.021 

8.15±0.29 

2.42 ±0.13 

23.94 ±0.90 

3.81±0.12 

0.800 ±0.003 1 



0.593 

8.15 

2.42 

24.68 

3.79 

0.799 



0.713 

3.39 

1.54 

11.52 

4.10 

0.787 


0.5 

1 0.716±0.046 

3.42 ±0.24 

1.49 ±0.14 

11.54±0.61 

4.01 ±0.28 

0.792 ±0.007 1 



0.733 

3.55 

1.42 

11.92 

4.13 

0.795 

7 


0.668 

6.46 

1.70 

22.07 

4.04 

0.789 


1 

1 0.679 ±0.013 

6.62±0.16 

1.68 ±0.06 

22.66 ±0.60 

4.14±0.11 

0.788 ±0.005 1 



0.680 

6.69 

1.66 

22.94 

4.14 

0.790 



0.581 

8.59 

2.06 

29.84 

3.46 

0.785 


1.5 

1 0.577 ±0.030 

8.72 ±0.43 

2.12 ±0.29 

30.69 ±0.98 

3.50±0.13 

0.781 ±0.014 1 



0.607 

8.92 

1.83 

31.66 

3.62 

0.795 



0.589 

2.95 

1.90 

10.48 

3.29 

0.793 


0.5 

1 0.593 ±0.029 

3.01±0.12 

1.85 ±0.07 

10.64 ±0.32 

3.30±0.19 

0.795 ±0.002 1 



0.568 

2.92 

1.89 

10.44 

3.13 

0.794 

8 


0.626 

6.42 

1.93 

22.24 

3.82 

0.794 

1 

1 0.673 ±0.047 

7.03 ±0.61 

1.81±0.12 

23.61 ±1.37 

4.19±0.37 

0.795 ±0.001 1 



0.669 

7.06 

1.79 

23.79 

4.14 

0.795 



0.709 

10.88 

1.71 

36.16 

4.64 

0.803 


1.5 

1 0.696±0.018 

10.83 ±0.30 

1.71 ±0.05 

36.22 ±0.70 

4.48 ±0.19 

0.803 ±0.002 1 



0.689 

10.84 

1.73 

36.38 

4.44 

0.802 


Table 9: Double-NBD fit parameters for charged-particle multiplicity distributions of INEL events. The shaded 
(middle) row for each entry contains central values of the parameters with their estimated uncertainty, and the 
surrounding rows represent the parameter values of the bounding fits: topmost in the first row and the bottommost 
in the last one. 


The shape evolution is quantified by the parameter ( 77 ) 2 , which tends to increase with increasing 77 range 
and with increasing centre-of-mass energy. The observed relation ( 77)2 w 3 x ( 77 ) 1 , is consistent with the 
analysis of CMS data reported in fTTl . despite the fact that the scaling parameter A was not used in their 
fit function. 


In the bins I 77 I < I and I 77 I < 1.5, the CMS NSD data = 0.9 TeV (Fig. 17 1 showed a different trend 
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VS 

(TeV) 

Ari 

a 

<n>i 

ki 

<n>2 

k2 

A 



0.317 

1.90 

4.47 

4.80 

2.45 

0.935 


0.5 

1 0.363 ±0.094 

2.07 ±0.23 

4.34±1.22 

5.31 ±0.65 

2.75 ±0.43 

0.930 ±0.014 1 



0.457 

2.30 

3.46 

5.92 

3.09 

0.922 

0.9 


0.548 

4.58 

2.80 

11.70 

3.62 

0.952 

1 

1 0.572 ±0.074 

4.98 ±0.57 

2.82 ±0.29 

12.79±1.41 

3.88 ±0.50 

0.944 ±0.012 1 



0.646 

5.55 

2.77 

14.19 

4.37 

0.937 



0.661 

8.28 

2.73 

19.28 

4.28 

0.933 


1.5 

1 0.703 ±0.120 

9.20±1.23 

2.54 ±0.24 

21.51±3.38 

4.73 ± 1.06 

0.937 ±0.012 1 



0.823 

10.43 

2.30 

24.89 

5.67 

0.936 



0.513 

2.49 

2.73 

7.42 

3.11 

0.934 


0.5 

1 0.509 ±0.019 

2.50 ±0.05 

2.77 ±0.28 

7.53±0.13 

3.04 ±0.09 

0.931 ±0.006 1 



0.493 

2.48 

2.66 

7.45 

2.95 

0.931 

2.76 


0.605 

5.64 

2.48 

16.36 

3.87 

0.941 

1 

1 0.545 ±0.060 

5.32 ±0.34 

2.71 ±0.24 

15.65 ±0.75 

3.52±0.35 

0.935 ±0.006 1 



0.498 

5.11 

2.87 

15.26 

3.24 

0.930 



0.703 

10.31 

2.38 

26.50 

4.63 

0.908 


1.5 

1 0.703 ±0.067 

10.63 ±0.92 

2.18±0.21 

27.39 ±1.77 

4.52 ±0.49 

0.913 ±0.006 1 



0.655 

10.04 

2.20 

26.47 

4.19 

0.915 



0.748 

3.87 

1.66 

12.24 

4.54 

0.944 


0.5 

1 0.688 ±0.060 

3.58±0.32 

I.81±0.15 

11.54±0.80 

4.07 ± 0.46 

0.942 ±0.002 1 



0.649 

3.38 

1.90 

11.16 

3.75 

0.940 

7 


0.673 

7.00 

2.02 

22.77 

4.29 

0.948 


1 

1 0.671 ±0.020 

7.07 ±0.27 

2.02 ±0.11 

23.06 ±0.74 

4.27 ±0.18 

0.942 ±0.006 1 



0.679 

7.25 

1.95 

23.49 

4.34 

0.943 



0.615 

10.17 

1.94 

31.42 

3.75 

0.931 


1.5 

1 0.622 ±0.019 

10.43 ±0.43 

1.89 ±0.07 

32.31 ±0.89 

3.84±0.13 

0.930 ±0.004 1 



0.618 

10.36 

1.89 

32.76 

3.82 

0.932 



0.557 

2.95 

2.10 

10.28 

3.18 

0.927 


0.5 

1 0.567 ±0.018 

3.03 ±0.08 

2.04 ±0.06 

10.48 ±0.20 

3.22±0.12 

0.929 ±0.002 1 



0.555 

3.01 

2.02 

10.44 

3.13 

0.929 

8 


0.589 

6.25 

2.25 

21.71 

3.68 

0.932 

1 

1 0.625 ±0.035 

6.70 ±0.44 

2.13±0.12 

22.72 ±1.00 

3.94 ±0.26 

0.932 ±0.001 1 



0.631 

6.88 

2.07 

23.18 

3.98 

0.931 



0.817 

14.18 

1.56 

40.81 

5.60 

0.917 


1.5 

1 0.857 ±0.042 

15.31± 1.14 

1.49 ±0.07 

43.06 ±2.25 

6.25 ±0.73 

0.918 ±0.003 1 



0.839 

15.18 

1.51 

42.89 

5.83 

0.915 


Table 10: Double-NBD fit parameters for charged-particle multiplicity disttibutions of NSD events. The shaded 
(middle) row for each entry contains central values of the parameters with their estimated uncertainty, and the 
surrounding rows represent the parameter values of the bounding fits: topmost in the first row and the bottommost 
in the last one. 

at high multiplicity as compared to the ALICE data (Section [9^ . We have made our own fit of CMS 
data excluding the bin Neb - 0 (Table [T^ and using the overall scaling parameter (d). As we cannot 
take into account correlations within CMS data, we compare these fits with fits of our final distributions, 
ignoring the correlations. We find that the relative weight of the second NBD component in the CMS 
data is smaller than in the ALICE data, while other parameters are compatible within their respective 
uncertainties. This confirms the different trends observed in distribution tails. 

In conclusion, a double-NBD function provides a precise description of the entire set of multiplicity 
distributions measured in this experiment. 
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9.6 KNO studies 

The KNO variable Ncul (Nch) provides another way to study the evolution of the shape of multiplicity 
distributions with varying centre-of-mass energies and varying pseudorapidity intervals. This study is 
carried out for the NSD event class only so that SD events, which may have a different behaviour, are not 
included in the data samples. 

The quantities of interest are derived from the original set of 72 multiplicity distributions and the resulting 
spread is used to estimate the various uncertainties presented in this section. 


9.6.1 Evolution of the shape of multiplicity distributions with Vs 


The KNO test in the range 0.9 to 8 TeV is limited by the multiplicity reach at 0.9 TeV. KNO-scaled 
distributions and their ratios were obtained for each of the available combinations of corrections with 
the same procedure used for multiplicity distribution measurements (averaging the 72 cases listed in 


Section 8.2 and using the spread as a measure of the systematic uncertainty). Bin to bin correlations 
were ignored when comparing KNO distributions and ^-moments at various centre-of-mass energies. 
Consequently, the relative errors obtained on the ratios are somewhat overestimated. Ratios between 
the two highest energies and 0.9 TeV exceed the value 2 at NaJ{N^h) larger than 5.5, 5 and 4.5, for 
\r]\ < 0.5, \r]\ < 1 and I 77 I < 1.5, respectively (Fig. [T9]). This confirms that KNO scaling violation increases 
with increasing pseudorapidity intervals. The shape of the KNO scaling violation reflects the fact that 
the high-multiplicity tail of the distribution increases faster with increasing energy and with increasing 
pseudorapidity interval than the low (Vch < 20) multiplicity part as already noted in Section |93| 


9.6.2 Quantitative KNO study with normalized q-moments 

Multiplicity distributions may be characterized by their normalized ^-moments defined as 


r = ^ 
{n)‘i 


(17) 


where ^ is a positive integer studied here for values 2, 3, 4 and 5, for NSD events (Fig.|20|and Table [T3]). 
For the three pseudorapidity intervals studied here (|? 7 | < 0.5, 1 and 1.5), C 2 remains constant over the 
energy range, C 3 shows a small increase with increasing energy for the two largest q intervals, C 4 and 
C 5 show an increase with increasing energy, which becomes stronger for larger q intervals. These new 



a 

<n>i 

ki 

<n>2 

k2 

A 

(TeV) 








0.595 

4.40 

2.56 

11.82 

3.69 

1.047 

0.9 [ 

0.661 ±0.069 

4.97 ±0.58 

2.45 ±0.20 

13.39±1.57 

4.21 ±0.55 

1.046 ±0.004 1 


0.731 

5.55 

2.25 

14.96 

4.76 

1.048 


0.581 

5.12 

2.37 

15.64 

3.57 

1.038 

2.76 [ 

0.588 ±0.028 

5.15±0.29 

2.30 ±0.11 

16.12±0.69 

3.59 ±0.24 

1.047 ±0.012 1 


0.572 

4.98 

2.19 

15.79 

3.42 

1.060 


0.671 

6.63 

1.72 

22.36 

4.13 

1.049 

7 1 

0.683 ±0.012 

6.65 ±0.24 

1.68 ±0.05 

22.86 ±0.55 

4.19±0.12 

1.053 ±0.010 1 


0.684 

6.78 

1.69 

23.08 

4.19 

1.052 


0.631 

6.26 

1.74 

22.07 

3.75 

1.056 

8 1 

0.640 ±0.009 

6.44 ±0.18 

1.73 ±0.04 

22.59±0.51 

3.84 ±0.09 

1.053 ±0.004 1 


0.641 

6.57 

1.76 

22.86 

3.85 

1.049 


Table 11: Double-NBD fit parameters for charged-particle multiplicity distributions of INEL>0 events. The 
shaded (middle) row for each entry contains centtal values of the parameters with their estimated uncertainty, 
and the surrounding rows represent the parameter values of the bounding fits: topmost in the first row and the 
bottommost in the last one. 
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VS 

Ai? 


a 

<W>1 


<w>2 

ki 

X^l ”dof 

(TeV) 








ALICE 

0.9 

1 

0.95 ±0.02 

0.56±0.10 

5.0±0.8 

2.6 ±0.8 

12±2 

3.4±1.1 

1.5/53 

1.5 

0.96 ±0.02 

0.70±0.18 

8.4±3.0 

2.5 ±1.1 

21±8 

4.7±3.1 

14.0/65 

CMS 

0.9 

1 

0.93 ±0.03 

0.72 ±0.08 

5.4±1.0 

2.6 ±0.1 

14±2 

5.8±2.3 

2.3 / 33 

1.5 

0.94 ±0.02 

0.74 ±0.05 

8.4±1.0 

2.6 ±0.1 

21±2 

6.4±1.8 

1. 8/45 


Table 12: Comparison of double-NBD fit parameters for charged-particle multiplicity distributions of NSD events 
by ALICE and CMS, obtained without accounting for correlations using the function defined by Eq. ( [T6| . CMS 
data are taken from fM). The last column gives the;^^^ value and the number of degrees of freedom (wdof)- 


data are in agreement with UA5 ETIl . CMS 1(3^ and ALICE’S previous results Elll in all cases where a 


comparison was possible (Eig. 201. 


10 Discussion of results and conclusion 

The AEICE Collaboration has carried out a detailed study of pseudorapidity densities and multiplicity 
distributions of primary charged particles produced in proton-proton collisions, at = 0.9, 2.36, 2.76, 
7 and 8 TeV, in the pseudorapidity range \t]\ < 2. A large increase of event sample size compared to 
previous AEICE publications, combined with improved measurement techniques, was used to study the 
evolution of charged-particle multiplicities over the whole centre-of-mass energy range covered so far 
by the EHC. The data at the highest energy appear as a smooth continuation of lower energy data, both 
in shape and in magnitude. 

The pseudorapidity density of charged particles, dAch/d?/, was measured as a function of pseudorapidity 


Vs (TeV) 


At] 

Moment 

0.9 

2.76 

7 

8 


(Ach> 

3.8±0.1 

4.6+ 0.2 

5.7+ 0.2 

5.8+ 0.4 


C2 

2.0±0.1 

2.0 + 0.1 

2.1+0.1 

2.1+ 0.2 

0.5 

C3 

5.3+ 0.4 

5.8 + 0.7 

6.1+0.5 

6.4+1.0 


C4 

19 + 2 

21+4 

23 + 3 

24 + 5 


C5 

78 + 11 

91+23 

102+17 

112 + 33 


(Aeh> 

7.8 + 0.3 

9.4+ 0.4 

11.6 + 0.4 

11.9 + 0.7 


C2 

L8 + 0.1 

L9 + 0.1 

L9 + 0.1 

2.0 + 0.1 

1 

C3 

4.3+ 0.3 

4.9+ 0.6 

5.3+ 0.4 

5.6+ 0.8 


C4 

13 + 1 

16 + 3 

18 + 2 

20 + 4 


C5 

48 + 7 

62 + 16 

73+12 

81+24 


(Ach> 

11.8 + 0.4 

14.2 + 0.7 

17.5 + 0.6 

17.8 + 1.1 


C2 

L7 + 0.1 

L8 + 0.1 

L9 + 0.1 

L9 + 0.1 

1.5 

C3 

3.9 + 0.3 

4.5+ 0.6 

5.0+ 0.4 

5.2+ 0.8 


C4 

11 + 1 

14 + 3 

16 + 2 

18 + 4 


C5 

36 + 6 

51 + 13 

64+ 11 

70 + 20 


Table 13: Mean charged-particle multiplicity and ^-moments of the multiplicity distributions measured by ALICE 
- 0.9, 2.76, 7 and 8 TeV, for NSD events in three different pseudorapidity intervals. The uncertainties are 
combined statistical and systematic uncertainties. 
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Fig. 19: KNO-scaled distribution (A^ch) P (iVch) versus the KNO variable A^ch/(iVch) at= 0.9, 2.76, 7 and 8 TeV, 
for three pseudorapidity intervals: \ri\ < 0.5 (top), 1.0 (middle) and 1.5 (bottom). In each case, ratios to the 
distribution at V* = 0.9 TeV are shown, on the right-hand side parts of the hgures. As Vch/(A^ch) takes different 
values at different centre-of-mass energies, ratios were obtained by interpolating the KNO-scaled distributions, and 
uncertainties were taken from the nearest data point. Bands represent the total uncertainties. 


in the range I 77 I < 2. The relative precision achieved at 77 = 0 for = 7 TeV is 5.5%, 2.6% and 1.3% for 
INEL, NSD and INEE>0 event classes, respectively. 

The power law parameterization of dAch/d ?7 at rj = 0, s^, provides a good description of the data from 
ISR to EHC energies: 6 = 0.102 ±0.003, 0.114 ±0.003 and 0.114 ±0.001, for the INEE, NSD and 
INEE>0 event classes, respectively, to be compared to d - 0.15 for Pb-Pb collisions |[54l . The AEICE 
Collaboration has shown clearly that the particle pseudorapidity density increases faster with energy in 
Pb-Pb collisions than in pp collisions. The extrapolation of dAch/d?; at 77 = 0 to the nominal EHC energy 
{^fs - 14 TeV) is obtained with a precision of 4.6%, 3.0% and 1.3% for INEE, NSD and INEE>0 event 
classes, respectively. 

Multiplicity distributions of primary charged particles were measured in three pseudorapidity ranges: 
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Fig. 20: Centre-of-mass energy dependence of the ^'-moments (g' = 2 to 4, left-hand scale, and q- 5, right-hand 
scale) of the multiplicity distributions for NSD events in three different pseudorapidity intervals (I 77 I <1.5 top, 
I;;! <1.0 middle and |? 7 | < 0.5 bottom). ALICE data (black) are compared to UA5 ETll (red) for \q\ < 0.5 and \q\ < 1, 
at V* = 0.9 TeV, and with CMS IlMl (blue) dXyfs — 0.9 and 7 TeV for \q\ < 0.5. The error bars represent the 


combined statistical and systematic uncertainties. The data at 0.9 and 7 TeV are slightly displaced horizontally for 
visibility. 


\j]\ < 0.5, \j]\ <1.0 and I 77 I <1.5. At = 7 TeV, N^h reaches about 70, 120 and 150, in I 77 I < 0.5, 1.0 and 
1.5, respectively, with the present statistics. These correspond to multiplicity densities 8 to 10 times the 
average multiplicity density. Based on the J. D. Bjorken formula |[83l . a characteristic collision energy 
density can be estimated, which increases by the same factor. For a qualitative estimate, assuming that 
the average energy density in pp collisions at = 7 TeV is of the order of 1 GeV/fm^ (see for example 
f/hl ). a density of 10 GeV/fm^ should be reached with high multiplicity pp collisions, similar to the 
energy density of Au-Au central collisions at RHIC If84l . When LHC runs at its nominal centre-of-mass 
energy of 14 TeV, high multiplicity proton-proton collisions will provide further direct comparisons of 
nuclear matter properties for interacting systems with similar energy densities but very different volumes. 

At LHC energies, above 2 TeV, multiplicity distributions can no longer be represented by a single NBD, 
but a double NBD gives a good representation of the data. The deviation from single NBD is already 
visible in the tail of the distributions = 0.9 TeV, but becomes increasingly large as the centre-of-mass 
energy increases. 

A test of KNO scaling between ^[s = 0.9 and 8 TeV confirms that KNO scaling violation increases with 
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increasing and, at a given centre-of-mass energy, with increasing width of pseudorapidity intervals. 

Comparisons with models in the pseudorapidity range I 77 I < 1 show that none of the event generators con¬ 
sidered is able to describe all properties of charged particle production up to = 8 TeV. PYTHIA 6 Peru¬ 
gia 2011, PYTHIAS 4C and EPOS EHC describe the pseudorapidity densities fairly well atV^ = 7 TeV, 
as well as the multiplicity distributions, but not above Vch ~ 60. The fact that PYTHIA 6 Perugia 2011, 
PYTHIAS 4C and EPOS EHC are in reasonable agreement with the data presented in this publication 
can probably be attributed to the fact that these generators were adjusted using the first EHC data. 

Studies of charged-particle production are refining the understanding of global properties of proton- 
proton collisions at the EHC. It was already demonstrated that, as with lower energy data, there is a 
strong correlation between multiplicity and mean transverse momentum (px) lES- However, there is also 
evidence 14] that the high multiplicity events have a topology more spherical than expected from current 
event generators PYTHIA 6 , PYTHIAS and PHOJET, suggesting that single hard-parton collisions may 
not be the only contributors to high multiplicity events. The general picture that emerges from this study 
is that, from ^/s = 0.9 to 8 TeV, multiplicity distributions and charged particle pseudorapidity densities 
follow a smooth evolution. 
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